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Introduction 



This Ph.D. thesis pursues two goals: The study of the geometrical struc- 
ture of two-dimensional quantum gravity and in particular its fractal nature. 
To address these questions we review the continuum formalism of quantum 
gravity with special focus on the scaling properties of the theory. We dis- 
cuss several concepts of fractal dimensions which characterize the extrinsic 
and intrinsic geometry of quantum gravity. This work is partly based on 
work done in collaboration with Jan Ambj0rn, Dimitrij Boulatov, Jakob L. 
Nielsen and Yoshiyuki Watabiki [1]. 

The other goal is the discussion of the discretization of quantum gravity 
and to address the so called quantum failure of Regge calculus. We review 
dynamical triangulations and show that it agrees with the continuum theory 
in two dimensions. Then we discuss Regge calculus and prove that a contin- 
uum limit cannot be taken in a sensible way and that it does not reproduce 
continuum results. This work is partly based on work done in collaboration 
with Jan Ambj0rn, Jakob L. Nielsen and George Savvidy [2]. 

In chapter 1 we introduce the main ingredients for the formulation of 
two-dimensional quantum gravity as an Euclidean functional integral over 
geometries. It contains a brief reminder of Liouville theory and the technical 
issues in the continuum formalism. We use these techniques to discuss the 
extrinsic and intrinsic Hausdorff dimension and the spectral dimension of 
two-dimensional quantum gravity. 

Chapter 2 is a review of dynamical triangulation in two dimensions. 
We begin with an introduction of the main ideas of how to discretize two- 
dimensional quantum geometries. The scaling properties are illustrated by 
means of the two-point function of pure gravity and of branched polymers. 

In chapter 3 we discuss quantum Regge calculus which has been suggested 
as an alternative method to discretize quantum geometries. We prove by a 
simple scaling argument that a sensible continuum limit of this theory cannot 
be taken and that it disagrees with continuum results. 
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Chapter 1 

Two-dimensional quantum gravity 



Euclidean quantum gravity is an attempt to quantize general relativity based 
on Feynman's functional integral and on the Einstein-Hilbert action princi- 
ple. One integrates over all Riemannian metrics on a d-dimensional manifold 
M. It is based on the hope that one can recover the Lorentzian signature af- 
ter performing the integration analogously to the Wick rotation in Euclidean 
quantum field theory. For a general discussion of further problems and for 
motivation of a theory of quantum gravity we refer to [3, 4, 5]. 

General relativity is a reparametrization invariant theory which can be 
formulated with no reference to coordinates at all. This diffeomorphism 
invariance is a central issue in the quantum theory. Its importance is most 
apparent in two dimensions, since the Einstein-Hilbert action is trivial and 
consists only of a topological term and the cosmological constant coupled 
to the volume of the spacetime. All the non-trivial dynamics of the two- 
dimensional theory of quantum gravity thus come from gauge fixing the 
diffeomorphisms while keeping the geometry exactly fixed. This is the famous 
representation of the functional integral over geometries as a Liouville field 
theory by Polyakov [6] . Based on this formulation the scaling exponents can 
be obtained [7, 8, 9]. 

Any theory of quantum gravity must aim at answering questions about 
the geometrical structure of quantum spacetime. The interplay between 
matter and geometry is well known from general relativity. The quantum 
average over all geometries changes the dynamics of this interaction. It turns 
out that the quantum spacetime has a fractal nature and even its dimension 
is a dynamical quantity. The characterization of the fractal nature of the 
quantum spacetime in two-dimensional quantum gravity is one of the central 
themes of this work. 

In section 1.1 the main concepts in the continuum formalism are intro- 
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duced. In the following section we review technical details about the factor- 
ization of the diffeomorphisms from the functional integral. In section 1.3 we 
introduce fractal dimensions to characterize the fractal nature of quantum 
spacetime. That section, and in particular section 1.3.2 are partly based on 
work presented in [1, 10]. 

Reviews about the continuum approach to two-dimensional quantum 
gravity can be found in [11] and the references therein. For an introduc- 
tion to conformal field theory see [12, 13, 14]. 

1.1 Scaling in quantum gravity 

In this section we introduce the basic concepts in two-dimensional quantum 
gravity. We begin with the partition function and the Hartle-Hawking wave- 
functionals. The most natural object to address questions about the scaling 
behaviour of the theory is the geodesic two-point function. We discuss its 
scaling behaviour in detail and introduce the concept of the intrinsic fractal 
dimension which illustrates the fractal nature of the intrinsic geometry of 
the two-dimensional quantum space-time. 

1.1.1 Partition function 

Let M be a two-dimensional, closed, compact, connected, orientable manifold 
of genus h. Then the partition function for two-dimensional quantum gravity 
can formally be written as the functional integral 

Z(G,A)=Y_ 

h>0 

Here 



d 2 £Vgft(£,) (1.2) 

M 

is the classical reparametrization invariant Einstein-Hilbert action [15, 16] 
with the gravitational coupling constant G, the cosmological constant A and 
the curvature scalar 1Z. According to the Gauss-Bonnet theorem, the last 
term in (1.2) is a topological invariant, called the Euler characteristic x of 
M: 



-S BH (g;G,A) 



VgX e - Smatter(9 ' X) . 



(1.1) 



s BH (g;G,A) = a 



d 2 ^n(i)=2-2\ 



(1.3) 
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while the first term is the volume V g of M equipped with the metric g^: 

d 2 V9. (1-4) 



M 



Therefore, (1.1) can be rewritten as 

X(K) 

h>0 

where Z(A) is defined as 



Z(G,A) = ^e^Z(A), (1.5) 



Z(A) 



V[g^] e - s{gA] V g X e - s ^uAg,x) t ^ ^ 



with S(g, A) = AV g . 

S m atter(g>X) in (l.l) denotes any reparametrization invariant action for 
conformal matter fields X with central charge D coupled to gravity. A typical 
example is the coupling of D free scalar matter fields X 1 , . . . , X D to gravity. 
In this case 



Smatter(g) X) — — 



d 2 Vg g^x^x", (i.7) 

M 



which is diffeomorphism invariant and invariant under Weyl rescalings of the 
metric: 

S m atter(e*g,X) =S matter (g,X). (1.8) 

Note that (1.7) can also be interpreted as an embedding of M in a D- 
dimensional Euclidean space, thus leading to an interpretation of (1.1) as 
bosonic string theory in D dimensions [6]. 

The functional integration J"D[g^] in (1.1) is an integration over all ge- 
ometries, that means all diffeomorphism classes [g HV ] of metrics g^ on the 
manifold M. This is often denoted formally as 



Vol(Diff) ' 



(1.9) 



where Vol(Diff) is the volume of the group of diffeomorphisms on M. Since 
this group is not compact the quotient (1.9) does not make sense beyond a 
formal level. We will derive expressions for the measures Vg^, Vig^] and 



V g X in section 1.2. 



In (1.1) the sum goes over all topologies of two-dimensional manifolds, 
that means over all genera h. It is presently unknown how to define such a 
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sum over topologies different from a perturbative expansion in h. It comes 
as a bitter aspect of the theory that we only know how to perform the 
functional integration in (1.1) for a given manifold and thus only for a given 
fixed topology, while a summation over topologies has to be performed by 
hand. In four dimensions it is known that the topologies even cannot be 
classified. Therefore the role of topological fluctuations is still unclear in two 
and in higher dimensions. In the remainder of this work we will always fix 
the topology and thereby disregard the sum over the genera in (1.1). 

Sometimes it is useful to define the partition function in an ensemble of 
universes with fixed volume V: 



Z(V) 



V[g^] 6(V-V g 



VgX e - Smatt " (9 ' X) 



VgX e - Smatt - (9 ' X) . 



Z(A) is the Laplace transform of Z(V) 1 , that means: 



Z(A) = 
Z(V) = 



dV e- AV Z(V), 



o 

rcr+ioo 



dA 
2m 



,AV 



Z(A), 



(1.10) 

(1.11) 
(1.12) 



where a is a real constant which exceeds the real part of all the singular- 
ities of Z(A). Formula (1.10) shows that although the action of pure two- 
dimensional quantum gravity is trivial and does not contain any propagating 
degrees of freedom, quantum gravity in two dimensions is a true quantum 
problem, since each equivalence class of metrics is counted once with the 
same weight. That means that there is no classical spacetime around which 
we expand and the theory is as "quantum-like" as it can get. Therefore 
it might contain important information about quantum gravity in higher 
dimensions. 

As we show in section 1.2 Z(V) scales as 



Z(V) oc V y -\ or Z(A) oc A 2 ^, 



(1.13) 



which defines a scaling exponent y that depends on the matter coupled to 
gravity and on the topology of the universe. The computation of y gives the 
result [7, 8, 9] 

1 -H 



y = 2 + — (D-25- V(25-D)(l-D)) 



(1.14) 



1 That is also the reason why we denote both with the same symbol and distinguish 
between them by the names of their arguments. 
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This formula is often called KPZ-formula. For pure gravity (D = 0) y equals 
—1/2, for the Ising model coupled to gravity (D = 1/2) we have y = — 1/3 
and at D = 1 , y equals for spherical topology. Clearly this formula breaks 
down for conformal charges D > 1 , where y assumes complex values. This 
has sometimes been called the D = 1 barrier. There is now some evidence 
for the fact that two-dimensional quantum gravity coupled to matter with 
central charge D > 1 is in a branched polymer phase [17, 18] in which the 
surfaces collapse to tree-like objects. 

The expectation value of some observable 0{g,X) is defined as 

V g X e" Smatter(9 ' x) C(g,X), (1.15) 

or for fixed volume V as 

V g X e- Smatter(9 ' x) C>(g,X). (1.16) 



(0(g,X)) A = 



ZfA). 



,-AV a 



|XVJ 



(0(g,X)) 



v 



Z(V) 



IxvJV 



1.1.2 Hartle- Hawking wavefunctionals 

Let us in the remainder of this section ignore the coupling of possible matter 
fields to quantum gravity. The definitions below can easily be expanded to 
the general case. Furthermore we concentrate on spherical topology only. 
Again the generalization is straightforward. 

Typical observables in two-dimensional quantum gravity are loop am- 
plitudes, that means amplitudes for one-dimensional universes. Let M be 
topologically equivalent to a sphere with b holes. The induced gravity on 
the boundary is covered by the modified action 

b 

S(g, A, Zi, . . . , Z b ) = AV g + Y_ ZiL g ,i, (1.17) 

i=1 

where L 9)i denotes the length of the i'th boundary component in the metric 
g^-v. In this case the partition function is given by 



W(A,Z 1 ,...,Z b ) = 



V[9w] e- s(9AZl Zt) . (1.18) 



This can also be interpreted as the amplitude for b one-dimensional universes 
of arbitrary lengths. Since the lengths of the boundary components are 
invariant under diffeomorphisms it makes sense to fix them to prescribed 
values L-! , . . . , L b . This yields the definition of the Hartle-Hawking wave 
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functionals [19]: 

W(A,L h ... ,L b ) = 



V[ 9]£V ] e- s( 9- A) nMU-L g ,0. (1.19) 



i=1 



We note that (1.18) is the Laplace transform of (1.19), that means: 



W(A,Z 1) ... ,Z b ) 



Y[dUe-^ W(A,L b ... ,L b ). 



(1.20) 



In the ensemble of universes with fixed volume V the wave functionals are 
given by: 



W(V,Li,...,L b ) = 



V[ 9iXV ] 5(V-V g )n 5 ( L i-L 9)i ). (1.21) 



i=i 



The cosmological constant A and the volume V are conjugate variables, that 
means: 



W(A,L 1( ... ,L b ) = 



dV e- AV W(V,L,,...,L b ). 



(1.22) 



In the case b = one recovers the partition functions (1.6) and (1.10). 
1.1.3 The two-point function 

To define reparametrization invariant correlation functions let d g (£,,£,') be 
the geodesic distance between two points £, and £/ with respect to the metric 
g^. Then the invariant two-point function is defined as 



G(A,R) 



M 



V[ 9irv ] e- s ^A) d^v/gUJ 6(R - d B (£,, £,')). 

(1.23) 



M 



G(A, R) can also be interpreted as the partition function for universes with 
two marked points separated by a geodesic distance R. The integrated two- 
point function, called the susceptibility x has the following behaviour: 



X(A) = 



dR G(A, R) = 



9 2 Z(A) 

aA 2 



A~ Y . 



(1.24) 



Therefore y is often called the susceptibility exponent. Yet another charac- 
terization of the exponent y from the branching ratio of the two-dimensional 
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universes will be explained in chapter 2. For fixed volume of spacetime the 
definition of the two-point function is 



G(V,R) = 



M 



d^^gJT) d^vWiMR-dgU^')), (1-25) 



M 



which is related to formula (1.23) via 



G(A,R) 



'00 

AV 







dV e- AV G(V,R). (1.26) 



Similar to (1.23) a whole set of invariant correlation functions which depend 
on the geodesic distance can be defined by multiplying the local measures 
d 2 E 1 ^ l /g^X) and d 2 E,'^g[E,') by powers of the invariant curvature scalars 1Z{E,) 
and W). 

Note that by (1.23) or (1.25) the concept of geodesic distance becomes 
meaningful even in quantum gravity. The two-point function can be inter- 
preted as the partition function of the ensemble of universes which have two 
marked points separated by a geodesic distance R. Both, its short distance 
behaviour and its long distance behaviour reveal the fractal structure of the 
most important metrics that contribute to the functional integral (1.1). 

An important property of the two-point function (1.23) is the inequality 

G(A,Rt +R 2 ) > const x G(A,Ri) G(A,R 2 ), (1.27) 

which follows from a simple gluing argument which is explained in chapter 2. 
Up to a constant, (1.27) is equivalent to the subadditivity of — logG(A, R), 
from which the existence of the limit 

lim " l0gG(A ' R) = M(A) (1.28) 

follows from general arguments (see for example [20]). Furthermore one can 
deduce that M(A) is an increasing function of A and M(A) > for A > 0. 
At this stage one cannot prove that M(A) scales to zero as A goes to zero. 
However let us assume that this is the case and verify it later by an explicit 
calculation of the two-point function in section 2.3: 

M(A) = c A v , (1.29) 

with a dimensionless constant c. That means that for large R » M(A) _1 the 
two point function falls off exponentially with a subleading correction: 

G(A, R) ~ A?~y e- cAVR , for R » (1.30) 
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The power correction can be found by applying dimensional arguments to 
(1.24). We conclude that the average volume of two-dimensional universes 
with two marked points separated by a geodesic distance R is proportional 
to R if R is large enough: 



<V> R = 



31ogG(A,R) 



3A 



R>M(A)- 



A^R. 



(1.31) 



That means, for large R typical universes have the shape of long thin tubes. 
On the other hand, for R ~ A _v the exponential decay turns over into some 
power law and we get 



(V) R ~R-, for R~ A 



(1.32) 



by inserting into (1.31). Per definition, the exponent of R in this equation 
equals the (grand canonical) intrinsic Hausdorff dimension dn = \- 

The large R behaviour of G(V, R) can be computed from (1.30) by a saddle- 
point calculation. The result up to power corrections is: 



G(V,R) ~ e 



■ r R 

, for — »1, 



(1.33) 



where c = — -(-vc) 1 -^. 

Another concept of fractal dimension can be applied in the ensemble of 
surfaces with fixed volume V for small distances R. Let us define with 



9 J 



M 



d 2 £VgUJ d 2 £,V^6(R-d B (£;,0) (1-34) 



JM 



the average length of boundaries on a manifold M with the metric g^, 
which have a geodesic distance R from a marked point. Then l(g,R)dR is 
the average volume of a spherical shell of radius R and thickness dR. The 
quantum expectation value in the ensemble of universes with fixed volume 
is given by 



(l(9,R))v 
Now we conclude from 



1 



Z(V) J 



^[g^lvUg.R) 



G(V,R) 
VZ(V) ' 



dR (l(g,R)) v = V 



(1.35) 



(1.36) 



that 



dim[V] = dim[R] + dim[(l(g, R)) v ]. 



(1.37) 
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Thus (l(g,R)) v has the scaling behaviour 

(l{g t W)v = V^*(j£) t (1.38) 

where F(x) is a function which falls off exponentially for large x, see (1.33). 
The (canonical) intrinsic Hausdorff dimension dn is now denned by the scal- 
ing of (l(g, R)) v ~ R dh_1 for small R. To be precise we expand (1.38) around 
R = and get: 

(l(g,R)) v ~ V 1_vdh R dH_1 , forR<V^. (1.39) 

For smooth d-dimensional manifolds we have dn = dh = d. If (l(g,R)) v 
stays nonzero and finite for V — > oo we have 

vd h = 1 , that means d H = d h . (1-40) 

This requirement is called the smooth fractal condition. It means that the 
average circumference of circles with a small geodesic radius R does not 
depend on the global volume of the universe. However it is well known that 
the smooth fractal condition is not fulfilled for the model of multicritical 
branched polymers [21], compare section 2.1.4. Therefore its validity is non- 
trivial and should not be taken for granted uncritically. It turns out that 
in two-dimensional quantum gravity the smooth fractal condition is fulfilled. 
That means that there is only one fractal dimension for short and for long 
distances together. 

For pure gravity, in the case of spherical topology, the two-point function 
G(A, R) can be calculated exactly with a transfer-matrix method [22] or 
alternatively by a peeling method [23] . The result 

~, * ^ . xia coshcA 1/4 R 

G(A, R) oc A 3/ = — - (1.41) 

sinh 3 cAV4R *> ' 

fulfills the standard scaling relations, but with non-standard exponents. 
G(A,R) falls off as e ~ 2cA ' /4R , from which we read off "V = 1/4, so that the 
Hausdorff dimension of pure two-dimensional quantum gravity equals four. 

Another scaling exponent, though not independent of y and -v, can be 
defined by: 

G(A, R) ~ R 1 ^, for 1 < R < A~\ (1.42) 

r| is called the anomalous scaling dimension. By expanding (1.41) for small 
R we get 

G(A,R) = 1-1aR + 0(R 3 ), (1.43) 
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and thus rj = 4. This is a notable result, since in ordinary statistical systems 
r\ is always smaller than 2. 

The so-called Fisher scaling relation 

y = y(2-T\) (1.44) 

relates the exponents defined above. It can be derived by applying dimen- 
sional arguments to (1.24). 

1.2 Liouville theory: A brief reminder 

A central problem in the continuum approach to two-dimensional quantum 
gravity is posed by the diffeomorphism invariance of the theory. It is com- 
paratively easy [24] to derive an expression for the formal integration Vg^ 
since one can define a natural scalar product on the cotangent space to the 
space of all metrics which defines a volume form in the same way as in finite 
dimensional Riemannian geometry. 

However, since the measure and the action are diffeomorphism invariant, 
these expressions are ill-defined. Therefore a gauge fixing and the factor- 
ization of the diffeomorphisms from the measure is required. This has been 
performed in [6] where the functional integration over geometries is expressed 
as a Liouville field theory. This theory has been developed and explained in 
[25, 26]. The measure for the Liouville mode is very complicated. Therefore 
two-dimensional quantum gravity has strictly speaking not been solved yet 
in the continuum approach. 

However, the critical scaling exponents could be obtained in the light-cone 
gauge [7] and later in the conformal gauge by consistent scaling assumptions 
[8,9]. 

1.2.1 Functional measures 

Let R(M) be the space of all positive definite Riemannian metrics on a d- 
dimensional manifold M and let T g R(M) be its cotangent space at a point 

9w e R(M). 

To define the functional measure Vg^ one makes use of the fact, that the 
cotangent space T g R(M) can be naturally equipped with a diffeomorphism 
invariant scalar product (-,-) T . This is used to define a measure V g bg^ on 
T g R(M) which descends to a measure on R(M). The line element ds 2 for 



12 



CHAPTER 1. TWO-DIMENSIONAL QUANTUM GRAVITY 



fluctuations 5g^ e T g R(M) of the metric can be written as 

ds 2 = (5g,6g) T = 



& d E,y/g Sg^G^Hg^, (1.45) 

M 



with the DeWitt (super)metric [27] 

Up to an overall normalization constant, (1.45) is the only diffeomorphism 
invariant, ultralocal, that means dynamically neutral distance on T g R(M). 
The constant C in (1.46) takes the value —2 in canonical quantum gravity 
[27]. In our framework C cannot be computed. The choice of C determines 
the signature of the metric (■, -) T as can be seen by splitting the fluctuation 
bg^ into its trace part g^6c and its tracefree part Sh^. This decomposition 
is orthogonal with respect to (-,-) T . On the tracefree subspace G^-* 13 has 
eigenvalue 1 , while it has eigenvalue 1 + on the trace sector. Thus, the 
DeWitt metric is positive definite for C > — | and negative definite for 
C < — |. For C = — |, G^' a(3 is the orthogonal projection on the tracefree 
part of 6g^. The inverse of the DeWitt metric, which satisfies 

G^ a ,G^= ](^K + KK), (1-47) 



is given by 



1 (XV 



a|3 — ^ I 9na9-vp + 9|xp9va- -[ + Cd 9^9«P I > (l- 48 ) 



as can be verified by multiplication with (1.46). 
The line element (1.45) can be rewritten as 



ds 2 = 



d d £, 



d d £/ Y_ T_ Sg^n^fo £/)Sg aP (£/), (1.49) 



M • ,vl |x<-v a<p 



with a^- a P(£„ £,') = v/glll r^-°*(£,) S(£,-£/), and the x ^±12 ma trix 
r^' a|3 (£,) is defined as 

4G^ a P(£,), if ^i<^ and cc < (3, 
r^^{Q = { 2G^' a|3 (£,), if ^i<^ and cc = (3, or ^i = ^ and cc < |3, 
G^' ap (£J, if ^i = ^ and cc = (3. 

(1.50) 
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This induces a volume form Vbg^ on the cotangent space, which descends 
to a volume form Vg^ on R(M): 

Vbg^ = v/detQM.-v,«p (£>) £7j Yl d6g^(£J, (1.51) 

£,eM, \i<v 

Vg^ = y/det £l^[E, y V) ]J d 9^)- (1-52) 

We first note that due to the 5-function the determinant of 0^' a|3 (£„ £,') 
factorizes into a product over all spacetime points £, e M. The local deter- 
minant of the ^11 x ^±1) matrix v/g(£j r^- a|3 (£,) can be evaluated by 
variation with respect to 51ogg = — g^Sg^: 

6 logdet (y/g = Tr 6 log (r^> a|3 ) + Tr 5 log (^g l^' ap ) 

= G^ aP 6G^+^±il61ogg 

= l(d+1)(d-4) Slogg, (1.53) 

where it has been used that the inverse of r^- ^ equals restricted to 

all indices with \x < v and a < (3. Therefore it follows that 

det( v ^r^)=Kg iA±I ^ 1 , (1.54) 

where the constant k can be determined by specializing to g^ = 6^, which 
gives k = (l + ^) /2 ' 2 \ Thus the measures (1.51) and (1.52) equal 

Vbg^ = const x J} J} g(£,)^d6g^(£), (1.55) 

f,eM n<-v 

Vg^ = const x [jjjg^^dg^), (1.56) 

f,£M |a.<-v 

which we normalize such that 

V&g^ e -^ 6s ' 5g) T = 1. (1.57) 

Note that cr = is the only exponent of g(£,) in (1.55) or (1.56) such that 
these measures are diffeomorphism invariant. With the measure 

YlYlg^Vdbg^), (1.58) 

the Gaussian integral (1.57) is resolved by computing the determinant of 
g(£,)2- 2 T^> a P. This can be done along the lines of (1.53). The result is 

det (g(^- 2 T^) = const x g(Q '^'V^ 4 ' , (i. 59 ) 
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which is a diffeomorphism invariant constant if and only if cr = ^jjf . The 
same result can be derived more rigorously by using the BRS-symmetry 
associated with general coordinate transformations [28, 29]. 

Analogously, the measures for scalar and vector fields can be derived. 
The corresponding metrics (•, -) s and (•, -) v are defined as 



(4>,*> s = 



d d £,^g cp(£,)cb(£,), (1.60) 

i 

d d £,v^^^(£,)g^(£,). (1.61) 



The measures T> g § and turn out as: 

V 9 <\> = Ylg(Ql (1.62) 

£,eM 

d 

1.2.2 Factorization of the diffeomorphisms 

By construction, the functional measure and the action in (1.1) and (1.10) 
are invariant under diffeomorphisms. Therefore, the measure overcounts 
physically equivalent configurations related by the group of diffeomorphisms. 
To avoid this overcounting of gauge equivalent metrics, the diffeomorphisms 
have to be factored from the measure. In two dimensions, this can be done 
in a covariant way in the conformal gauge. 

Let us first recall that the space of metrics on M modulo diffeomorphisms 
and Weyl transformations is a finite dimensional compact space, called the 
moduli space M. of M, which has dimension 0, 2 and 6h— 6 for genus 0, 1 
and h > 2 respectively. It is parametrized by the Teichmiiller parameters 
t = (fi, . . . , t N ). That means that if for each moduli t e M. we choose a 
fixed background metric gV^M, all other metrics are contained in the 
orbits under diffeomorphisms and Weyl transformations: 

g^ = f*( e %v(T)), C 1 - 64 ) 

where f* denotes the action of the diffeomorphism f : M — > M. (Actually 
the background metrics 0^ can be chosen such that they all have constant 
curvature.) 

Any infinitesimal change 5g^ in the metric can be decomposed into an 
infinitesimal Weyl transformation S4>, the effect of an infinitesimal diffeo- 
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morphism £, and the effect of varying the Teichmiiller parameters: 

bg^ = f*6cj) + (V^ + V^) + ^^i- (1-65) 

Here denotes the covariant derivative. We want to orthogonalize this 
decomposition with respect to the scalar product (•, -) T . First note that the 
tracefree part of the effect of the diffeomorphism £, is given by the conformal 
Killing form 

(P£V-v = + V^-g^V a £, a , (1.66) 

which maps vector fields into traceless symmetric tensors. The zero modes 
of P are called conformal Killing vectors. For spherical topology there are 
six linearly independent conformal Killing vectors, for torus topology there 
are two while there are none for higher genus. The conformal Killing vectors 
are important because they induce a variation Sg^ in the direction of g^ 
which can thus be compensated by a Weyl rescaling. In other words, the 
decomposition (1.65) is in general not unique. To make it unique one chooses 
a diffeomorphism £, e (kerP)^ orthogonal to the zero modes of P. That 
means the gauge will only be fixed up to the conformal Killing vectors and 
we expect a remaining, at most six-dimensional symmetry in the gauge-fixed 
expressions. It follows that all variations 5g^ G T g R(M) of the metric can be 
uniquely decomposed into the variation of the conformal factor, the action 
of a diffeomorphism orthogonal to the conformal Killing vectors and the 
variation of the Teichmiiller parameters. This decomposition can be written 
as a mutually orthogonal sum of a trace part, a tracefree part orthogonal to 
the moduli deformations and these moduli deformations: 



+ (1-P(P+P)-ip+)lcJ rv 6T i , (1.67) 



6 9(ev = {f + 5c( ) + V cc l a +-g ap ^5T i }g^ + {P^ + P(P + P)- 1 P + k^T l } 



where 



The adjoint P+ of P is defined by the relation (h, P£,) T = (P+h, £,) v . To prove 
that the decomposition (1.67) is orthogonal just note that 1 — P(P + P)~ 1 P + is 
the projector on the zero modes ijn of P + . The tracefree part is orthogonal 
to the trace part by definition of the scalar product (•, -) T . 
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The change of variables g^ — > (4>,f,T) involves a Jacobian J(<J),t): 



V[gJ\ F(g) = 



Vf 



Vf 



(1.69) 



Here one integrates some reparametrization invariant functional T. To find 
J one substitutes the orthogonal decomposition (1.67) into the normalization 
condition 



1 



(1.70) 



and performs the change of variables on the tangent space. The result of 
these Gaussian integrations is 



1 =J(*,T) 



det (xhc,i|n)T 
det'(P+P) 



det ( t|> 



3g 



3t 



T 



-1 



(1.71) 



At last, the change from the diffeomorphisms f to f and the conformal Killing 
vectors a> a , which fulfill P(e*cu Q ) = 0, has to be computed. The final result 
is 



v[g^} Hg) = 



Vf, 



v e <t> Q 4) 



dT , / , 3g 
— - - det ( \p m , — ) x 
v(t) \ 3T n / T 

det'fP+Pl 



det (tp k) t|>i) T det (a) a , cu b ) 



v. 



■F(e*gtr)), (1.72) 



where v(t) is the volume of the group generated by the conformal Killing 
vectors. Now the diffeomorphisms can be factored out. Since it turns out 
that ^ det /"4> m , depends only on the Teichmiiller parameters and 

not on (|) we denote this as Vx. The square root is called the Faddeev-Popov 
determinant and denoted as Zppfe*^). 



1.2.3 Gravitational dressing of scaling exponents 

After factoring the diffeomorphisms from the measure the gauge-fixed ex- 
pression for the partition function (1.10) is: 



Z(V) 



V e ^<$> 6(V- V e 4. 9 ) Z F p(e*a)Z ma t(e*a), 



where Z mat is defined as: 



Z mat (e*g) 



(1.73) 



(1.74) 
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with some matter action for matter with central charge D. The measure for 
the Liouville mode 4> is defined by the line element 



d 2 s = (S4),S*>: 



d^v^encM^ScMa (1-75) 

M 



which depends on $ itself. Therefore the corresponding measure is 
very complicated and it is unknown how to perform the functional integration 
over § with this measure. The idea used to overcome these problems is to go 
over to the translational invariant measure V§ty and to shift all dependence 
on 4) in (1.73) into the action. 

For this transformation we use the relations 

Z mat (e*§) = Z mat (a)e^ SL( ^', (1.76) 
Z FP (e*3) = Z FP (3)e-^ s ^.<M (i.77) 
= e^^'* 5 , (1.78) 

where S L (§, cb) is the Liouville action: 

Sl(0, *) = \ m d2 tV§ Qr^M) + • (1-79) 

Here it is the curvature scalar in the metric 

The relation (1.76) for the matter partition function can be derived by 
a couple of methods such as the heat kernel regularization [30] or the zeta 
function regularization [31]. Since the matter action is invariant under Weyl 
transformations we can actually say that the whole contribution comes from 
the transformation of the measure gX under conformal rescalings of the 
metric. This is also the reason for the trace anomaly of the energy momentum 
tensor. 

The transformation property (1.77) for the Faddeev-Popov determinant 
can be derived in a similar way. However, the involved operators are much 
more complicated, compare (1.72). 

The transformation (1.78) could not be derived rigorously by these or 
other methods. David [8] and Distler and Kawai [9] use a selfconsistent 
bootstrap method to find this relation. Under the assumption that the fac- 
tor takes the form of an exponential of a Liouville action with arbitrary 
coefficients the coefficients are determined from general invariance consider- 
ations. 

In this calculation the Liouville mode § is rescaled by a factor <x, which 
has to be determined. Taking all factors together, the partition function 
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(1.73) can be rewritten as: 



Z(V) = 



Zpp(3) Z mat (§) e^ s ^>4» 6(V-V ea ^) ( (1.80) 



where the measures depend only on the background metric g\ Thus the 
partition function scales under a rescaling of the volume as: 



Z(AV) =A^ X - 1 Z(V), 



(1.81) 



where x is the Euler characteristic of the manifold M. To find (1.81) we have 
shifted 4) by ^logA and used the Gauss-Bonnet theorem (1.3). 

Now we have to compute the dressing exponent a. It turns out that 
more generally we can compute the anomalous scaling dimensions for a whole 
family of operators. Consider the observable 



0{g) 



dVg <D(g), 



M 



(1.82) 



where O is a spinless primary field of dimension A . Typical simple ex- 
amples are the volume V g = j M d 2 h,y/g with A = and the identity 1 = 
j M d 2 ^y/g -^=5(£,) with A = 1. Then (9(g) transforms under a rescaling 
g^-v — > Ag^ of the metric as 



(1.83) 



if the geometry is kept fixed. Under the average over all geometries the 
scaling behaviour exponent A is in general changed to some other value A: 



(O(g)) 



1-A 



AV 



(O(g)) 



v • 



(1.84) 



The coupling to gravity dresses the scaling behaviour of the field <D. There- 
fore we make the ansatz 



M 



d 2 £Vg e^O(g) 



(1.85) 



for the observable after the transition to the translational invariant measure. 
Note that 0(0, <J)) is diffeomorphism invariant and fulfills O{g,0) = 0{§). 
Now observe that (1.76), (1.77), (1.78) and 0(e*§) are invariant under the 
transformation 



(1.86) 
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Therefore also (9(§,(J)) has to be invariant. This requirement leads to an 
expression for the gravitational dressing exponent (3 as follows. Under (1.86) 
the product transforms as 



Thus e 13 *^' has to transform as 



e (Ao-1)ff(£.) e P()>(£,)_ 



(1.87) 



(1.88) 



On the other hand we can determine the transformation of e 13 *^' by com- 
puting the Gaussian integral 



D 9 4) e 



D-25 



(1.89) 



where |3<J)(£,) is interpreted as an additional 6-function source. If we define 



J(x) = H{x) + 



48tt(3 



1 



D-25 V /§M 



S(x-£,), 



(1.90) 



we get: 

F^(§,|3) = exp 



D-25 



96tt 



d 2 x^gtW(Ag 1 T)(x) 



M 



(1.91) 



M§,0) exp 



P 



M 



Under the transformation (1.86) the first term in the exponent gives an extra 
— |3cr(£,). The second term is an evaluation of the propagator at coinciding 
points which naively seen is singular. Thus we have to renormalize this 
expression [32] . The renormalized propagator is defined as 

Agite, Q = lim (Ag 1 (x, £,) - ^ log dg(x, ^ , (1.92) 

where dg(x, £,) is the geodesic distance between x and £,. Therefore the 
second term in the exponent gains an extra j^-^<j(Q under (1.86). Taking 
both contributions together we conclude: 



,(34>(x) 



e -(3o-(x)+^o-{ X ) e |3c|>(x) ± g(Ao-1)a(x) g |3c|>(x)_ (1.93) 

The solution to this quadratic equation determines (3 and thus also a: 

(3 = ^ (25 - D - y/(25 — D)(l — D + 24A )) , (1.94) 
a = ^(25-D-V(25-D)(l-D)). (1.95) 
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The sign of the square root is determined by requesting that in the classical 
limit c — > — oo the scaling of observables is the same as in a fixed geome- 
try. With <x inserted into (1.81) we can now verify the form (1.14) for the 
susceptibility exponent y. 

For the expectation values of observables the dependence on the topology 
cancels , and we get: 

(0(g)) xv = ^ (0(g)) y = A ] - A (0(g)) v . (1.96) 

A evaluates as 



A=V1 -U ±£ 4A -y^D 

V25 - D - yT^D v 1 

The expectation value of products of observables scales with the sum of their 
scaling exponents. 



1.3 Fractal dimensions 

Quantum gravity as defined above is a functional integration over geometries. 
Therefore we have to ask which geometrical concepts survive the quantum 
average and how geometry is changed by it. In any theory of quantum 
gravity the structure of spacetime is of primary interest and even its di- 
mensionality is a dynamical quantity. The transfer matrix method in [22] 
applied to pure two-dimensional quantum gravity revealed that the space- 
time has a self-similar structure on all scales with a fractal dimension four. 
An understanding of the quantum spacetime in two-dimensional quantum 
gravity might help to understand similar dynamical changes of the geometry 
in higher dimensions. 

In this section we introduce three different concepts of dimensionality 
to characterize the fractal structure of spacetime, each of which might de- 
scribe different aspects of the geometry. The intrinsic Hausdorff dimension 
or fractal dimension dn has already been defined above by the scaling of ob- 
servables such as the average circumference of circles with their "geodesic" 
radii. These measurements are performed within the intrinsic geometry of 
the universes. They are analogous to the usual "clock and rod" measure- 
ments in classical gravity. General results for d H from Liouville theory have 
been given in [33]. 

We can also characterize the fractal structure of spacetime by embedding 
it into a D-dimensional, typically Euclidean flat space. Measurements of, 
for example, the mean square extent of the universes are then performed 
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with the metric of the embedding space. The finite size scaling of the mean 
square extent of the surfaces defines the extrinsic Hausdorff dimension D H . 
Intuitively the embedding gives us an idea what typical universes would 
"look like". Indeed it is known that the extrinsic Hausdorff dimension of 
two-dimensional quantum gravity is infinity [34, 35]. That means that the 
surfaces are crumpled together in a complicated way. 

Instead of determining the metrical structure of the spacetime directly by 
putting out clocks and rods, one can also choose to observe the propagation of 
test particles. Prom such observations one can in the classical case determine 
the propagator and via an inverse Laplace transformation the associated heat 
kernel. Since all this can be done in a reparametrization invariant way, we 
can define corresponding quantum observables. This leads to the definition 
of the spectral dimension d s as the scaling exponent of the short time return 
probability distribution of diffusing matter. The first general analytical result 
for the spectral dimension of two-dimensional quantum gravity coupled to 
conformal matter with central charge D was given in [1]. 



1.3.1 Extrinsic Hausdorff dimension 



The extrinsic Hausdorff dimension D H of a fractal object embedded in a D- 
dimensional space is a measure for the extent of the object measured with 
the metric of the embedding space. D H can be defined for mathematical 
fractals as well as for random surfaces and for quantum gravity. Let us focus 
on quantum gravity. The partition function for two-dimensional quantum 
gravity with fixed volume V of spacetime coupled to D Gaussian matter 
fields is 



Z(V) 



IxvJV 



£> 9 XJ e 

y r-\ cm 



(1.98) 



Here V a X u \ 

» H cm 



denotes the functional integration over the D matter fields 
such that the centre of mass is fixed to zero, which amounts in dropping the 
zero modes of the Laplacian when integrating over the matter fields. This is 
a slight variation of (1.10) which allows us to define the extrinsic Hausdorff 
dimension as 



<X 2 ) V ~ V^, for oo, 



(1.99) 



with the definition 



1 



Z(V) 

1 



x 



DV 



^[g^V 

2r r= v2 



£>Xj e 

y H cm 



M 



d 2 £Vg X£(«, 



(1.100) 
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which is equivalent to 



<* 2 > 



1 



v 



Z(V) J 

1 



x 



2DV 2 



^[gjrvlv 



£> g Xj e" 

y H cm 



M 



d 2 £, 2v /g (X^-X^)) 2 . (1.101) 



M 



To understand the definition of D H note that the fields define an embed- 
ding of the manifold M in a D-dimensional Euclidean space. That means 
that we measure the Euclidean mean-square extent of spacetime embedded in 
D dimensions. Mathematically more conservatively the extrinsic Hausdorff 
dimension is defined by constructing a covering of the geometrical object 
embedded in some D-dimensional space by a union of n small balls, n de- 
pends on the typical macroscopic scale r of the system and scales with the 
dimension D H and is related to the mean square extent of the system as 



n 



(1.102) 



The extrinsic Hausdorff dimension for two-dimensional quantum gravity 
has been computed in [34, 35]. Note that the treatment in [34] is based on a 
more general scaling assumption than is needed in two-dimensional quantum 
gravity. Their general result was specialized to quantum gravity in [35]. In 
the following we briefly summarize the derivation. 

If we define the two-point function by 



G(p) = 
it follows that 



d 2 £, lv /g 



M 



a 2 Wg e 



ip{X(£,,)-X(£, 2 )) 



M 



V 



i a 2 

<X2 >v = -2DV^3p^ G(p) lp^ 



(1.103) 



(1.104) 



In flat space the two-point function behaves as V 2 - a "(p) with A ~ p 2 . Cou- 
pled to gravity this acquires a gravitational dressing which can be computed 
to be: 



A(p) 



y'l -D+24Ao(p)- y/T^D 
V25-D - ^/T^D 



(1.105) 



If we differentiate G(p) ~ V 2 A(p) two times with respect to p at p = we 

V— >oo 

get 



<x 2 > 
<x 2 > 



logV, for D < 1, 
log 2 V, for D = 1, 



(1.106) 
(1.107) 
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for large volumes V. In both cases (X 2 ) v grows slower than any power of 
V. Thus we conclude that D H = oo. For spherical topology this has already 
been found for D — > — oo and to one-loop order in [36] . 



1.3.2 Spectral dimension 

To probe the structure of classical spacetime one can study the properties 
of propagating particles, typically of free particles or diffusing matter. Let 
¥(£,,T) be the wave function for diffusion on the d-dimensional compact 
manifold M. W depends on the time-parameter T and is a function of the 
points in M which fulfills the diffusion equation 

■jfVfaV = A fl Y(fc,T) (1.108) 

for the initial condition 

V(t,0) = -J=6fo-y. (1.109) 

Here A g is the Laplace-Beltrami operator corresponding to the metric of 
M. The solution of the diffusion equation can be written as 



¥(£,,T) = e TA ^(£,,0) 



M 



d d tVgKi(t.t';TMt',o), (1.110) 



which defines the probability distribution (or heat kernel) K g (£,,£/;T) for 
diffusion on a compact manifold M with metric g^. Here and in the following 
we take into account, that the Laplace operator A g on compact surfaces has 
zero modes which should be projected out. This is indicated with a prime. 
Kg(£,, £/;T) is related to the massless scalar propagator (— A g ) _1 by 



(i'\(-A 9 r%y = 



dTK g (£,,£/;T). (1.111) 

o 



It is known that the average return probability distribution 

d<Vg K g (£,,£,; T) (1.112) 



M 



at time T for diffusion on a d-dimensional manifold with a smooth geometry 
[g^-v] admits the following asymptotic expansion for small times T: 

1 00 1 
Rp g(~n -jviY- a - jn = t^ (1 + ° (T]) ' (L113) 

r=0 
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where the coefficients a n are diffeomorphism invariant integrals of polyno- 
mials in the metric and the curvature and in their covariant derivatives, see 
for example [37] . This asymptotic expansion breaks down when T is of the 

2 

order V* , when the exponential decay of the heat kernel becomes dominant. 

Since RPg(T) is reparametrization invariant, we can define its quantum 
average over geometries for fixed volume V of spacetime as: 



RPv(T) = 3 



Z(V) J 



^[g^]ve- s ^ (s) RP:(T). (1.114) 



Here S eff (g) denotes the effective action of quantum gravity after integrating 
out possible matter fields. 

The spectral dimension d s is now defined by the short time asymptotic 
behaviour of RPy(T): 

RPvm-^n + om). (i.ii5) 

It is natural to expect that under the average (1.114) over all geometries the 
only remaining geometric invariant is the fixed volume V. Thus we expect 
that RPy(T) has the form 

RP v(T) = -^F (^) , (1.H6) 

where F(0) > and F(x) falls off exponentially for x — > oo. This scaling 
assumption is the main input into our derivation of the spectral dimension 
of two-dimensional quantum gravity. It is very well documented in numerical 
simulations [38, 39]. 

For a fixed smooth geometry the spectral dimension is by definition equal 
to the dimension d of the manifold. However, the quantum average can a 
priori change this behaviour. Actually we know that the intrinsic Hausdorff 
dimension of two-dimensional quantum gravity is different from two. In this 
sense generic geometries are fractal with probability one. For diffusion on 
fixed (often embedded) fractal structures it is well known, that the spectral 
dimension can be different from the embedding dimension as well as from 
the fractal dimension. The diffusion law, measured in the embedding space, 
becomes anomalous: 

<r 2 (T)) ~T 2AU , (1.117) 

with a gap exponent d w > 2. This slowing down of the transport is caused 
by the fractal ramification of the system. The gap exponent d w is related to 
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the spectral dimension d s and the intrinsic Hausdorff dimension dn by 

, 2d H 



(1.118) 



If the diffusion law is not anomalous, we have d w = 2 and d s = d h , analogous 
to d s = d for diffusion on fixed smooth d-dimensional geometries. For a 
review of diffusion on fractals see [40]. 

While the above concepts are in principal valid for Euclidean quantum 
gravity in any dimension d, let us now specialize to two dimensions. The 
starting point is again the partition function (1.98) for two-dimensional quan- 
tum gravity coupled to D Gaussian matter fields X,^ for fixed volume V of 
spacetime. We use (1.100) to define (X 2 ) y . The Gaussian action implies 
that: 



<x 2 > 



V 



DV3cu\ /v 
1 3 



a)=0 



DVZ(V) duo 
1 



£>[g^]v (det'(-A g - a))) 



-D/2 



tu=0 



2VZ(V) J 



V[g^] v (det'(-A g )) D/2 Tr' 
— 1^ • 

.-AgJ/y 



= —(TV 



-A, 



(1.119) 



By inserting (1.111) and (1.116) into this formula we get: 



(x 2 > 



V 



1 

2V 
1 

2 



dT 



d 2 ^ K' g (t, f,;T) 

M / V 



1 

dT RP^(T) = - 



1 T 



for V — > oo. Comparing this with (1.99) we conclude: 



J_ \ 

D^ + 2' 



Using the result D H = oo from the preceeding section we arrive at: 

d s =2, for all D < 1. 



(1.120) 



(1.121) 



(1.122) 



Strictly speaking, this assumes d s < 2. However, at D = — oo a fixed ge- 
ometry implies d s = 2 and D H = oo, in agreement with (1.121), and we 
expect that a saddle point calculation around D = — oo is reliable, compare 
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[36]. Therefore d s equals 2 in a neighbourhood of D = — oo. The average 
over all geometries in (1.114) includes many degenerate geometries, for which 
D H < oo, for example the branched polymer-like geometries which are dis- 
cussed in section 2.1.4. Thus, a priori we would expect that under the average 
over fluctuating geometries the spectral dimension decreases. Therefore it is 
reasonable to assume that d s < 2. 2 

Thus we have shown, that the spectral dimension of two-dimensional 
quantum gravity coupled to D Gaussian fields equals two for D < 1 [1]. Our 
derivation relies on the scaling form (1.116) of the averaged return probability 
RPy(T) and on properties of Gaussian matter fields. To corroborate the result 
(1.122), let us present another derivation completely within Liouville theory 3 . 
Let us define an observable 

1 1 



0( 9 ) = 



M 



-A fl yfilQ 



(1.123) 



Under a rescaling of the metric, O behaves like 0{\g) = AO(g), that means 
it is a correlator with conformal weight (—1 , — 1 ). Therefore we deduce from 
(1.96) that 



(0(g)> AV = M0(g)> 



v > 



or equivalently: 



On the other hand we can write: 



(1.124) 



(1.125) 



AV 



V 



0(g] 



v 




d 2 £, 0V / gW -\ 

M — A 



V 



d^ox/gl^I 

dT RPy(T) ~ 



dT e 




using the scaling assumption (1.116). From (1.125) and (1.126) it follows 
that d s = 2 in Liouville theory for all types of matter coupled to gravity 

2 If d s > 2, the integral in (1.120) diverges at 0, and we have to introduce a cut-off e 
for small times T. Then it is convenient to consider ^(X 2 ) n ) y ~ V 2ti ^ Dh instead, with 
n = [d s /2] + 1 for non-integer d s . Then the leading large V behaviour will be V 2n / d s^ 1 j 
and we get J- = ^ + ^ . 

3 Thanks to Jakob L. Nielsen for showing me this argument. 
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with central charge D < 1 , at which point the continuum calculations break 
down. 

This situation is very remarkable: A generic geometry in the functional 
integral (1.98) is a typical fractal, when looked at in the usual way by com- 
puting the fractal dimension d h from measurements of the circumference 
(l(g,R)) v of circles versus the geodesic distance R. Also the gap exponent 
d w is anomalous and larger than two. It is exactly equal to the intrinsic Haus- 
dorff dimension and consequently, the spectral dimension d s equals two. It 
takes the same value as for smooth two-dimensional geometries. 

The spectral dimension has received considerable attention also in nu- 
merical studies of the fractal structure of two-dimensional quantum gravity. 
In [39] the return probability is computed for pure gravity and for gravity 
coupled to matter with central charge \ and 1 . They find that their results 
are consistent with d s = 2. In [38] d s is measured for central charge —2, 0, \ 
and \. They get the results d s = 2.00(3), d s = 1.991(6), d s = \ and 
d s = 1.991(5), respectively 4 . 

For D > 1 it is generally believed that two-dimensional quantum grav- 
ity is in a branched polymer phase, see [17, 18] for recent analytical and 
[41] for some numerical evidence. The Gaussian fields define an embed- 
ding of these polymers in M D . The extrinsic Hausdorff dimension of generic 
branched polymers is D H = 4, thus we conclude from (1.121) that the spec- 
tral dimension equals |, the famous Alexander-Orbach value [42]. The value 
d s = | and formula (1.121) for branched polymers have been derived in [43] 
by different methods, see also [44] for a recent complete proof of d s = |. 
Furthermore, for the multicritical branched polymers [21] it is known that 
D H = -^j-, m = 2,3,..., where m = 2 for the ordinary branched poly- 
mers. Thus we obtain d s = 2 ^ a ; 1 , in agreement with the analysis in [45]. 
For m — > oo the multicritical branched polymers approach ordinary random 
walks, for which D H = 2 and d s = 1 . 

Let us in the end comment on a subtlety in the definition (1.114) of the 
averaged return probability RPy(T). We have defined it as the quantum aver- 
age of the return probability for fixed geometry over fluctuating geometries. 
Instead, we could have defined it as the limit for R — > of 



4 Thanks to Konstantinos N. Anagnostopoulos for showing me these data prior to the 
publication of [38] 



K V (R,T) = 




(1.127) 
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where G(V, R) = V (l(g, R)) V Z(V) is the partition function for universes with 
two marked points separated by a geodesic distance R ( "two-point function" ) 
as defined (up to matter fields) in (1.25), see also (1.35). Ky(R,T) is the av- 
erage probability distribution for diffusing a geodesic distance R in the time 
T. It is natural to identify RPy(T) with K v (0,T). However it is unknown 
whether the limit R — > of (1.127) commutes with the functional integration 
over geometries. Actually, for the two-point function G(V, R) it is easy to see 
that this is not the case. If the limits do not commute there are two inequiv- 
alent definitions of the return probability RP V (T). However, the scaling form 
of K v (0,T) is the same as the scaling form (1.116) of RP V (T) [46]. Indeed, 
we have 

'00 

dR(l(g,R)) v K V (R,T) = 1, (1.128) 

L 

and using dim[R] = dim[V v ] we arrive at the scaling form 

KWR,T) = lp(A, T.), (1 . 129) 

where A is chosen such that is dimensionless. Expanding the return 
probability Ky(0,T) around T = we get 

K v(0. T )~-fi7^> forT~0, (1.130) 

by definition of d s . If this short time asymptotic stays nonzero and finite for 
infinite volume V we get A = ^, and K v (0,T) scales as 

K v(0>T) = ^f(^), (1.131) 
which has the same form as the scaling behaviour of RPy(T). 

1.3.3 Intrinsic Hausdorff dimension 

The intrinsic Hausdorff dimension dh of two-dimensional quantum grav- 
ity has been defined in section 1.1.3. One way is to measure the average 
circumference (l(g,R)) v of circles with a small radius R, which scales as 
(U9)R))v ~ R dh ~ 1 - This is analogous to the classical "clock and rod" pro- 
cedure to study the geometrical properties of spacetime. Alternatively one 
can define the fractal dimension dn as the scaling exponent of the average 
volume of the universes: (V) R ~ R dH , compare (1.32). In two-dimensional 
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quantum gravity the first "local" and the second "global" concept of the 
intrinsic Hausdorff dimension agree. For pure gravity, we have d H = 4, see 
(1.41). 

The intrinsic Hausdorff dimension has been computed in [33] by study- 
ing the diffusion equation in Liouville theory. Although the mathematical 
status of this derivation is unclear and although it is not consistent with the 
derivation of the spectral dimension in section 1.3.2, the result is valid for a 
number of special cases and agrees with many numerical simulations. 

In [33] the mean square distance (R 2 (T)) y of diffusion in the fluctuating 
space time is defined as: 

1 



<R 2 m> 



V 



X>[q,Jv e~ Seff(9) 

VZ(V) J Ly ^ JV 



cn lv /g 



M 



M 



(R 2 (T)) V scales with the volume as 

(R 2 m 



v 



(1.132) 



(1.133) 



which yields yet another way to define the intrinsic Hausdorff dimension d H . 
On the other hand, the mean square distance for diffusion on a fixed smooth 
geometry [g^] with the initial condition (1.109) can be expanded for small 
times T as [47] 



M 



d 2 £V9 d 2 Jt , £,)K'(£,o, t;T) = 4T - -l 2 n{l ) + 0(T 3 ), (1.134) 



with the curvature scalar 1Z. If we assume that this expansion commutes 
with the functional average over geometries in (1.132), we get 



(R 2 (T)) V ~ const • T, for small times T. 



(1.135) 



The dimension of T can be computed from Liouville theory by expanding 
the (unnormalized) return probability for small times to observe how this 
scales under the transformation V — > AV: 



d 2 Eoy/g Ki(£, ,£,o;T)) = 



M 



AV 



M 



M 



d 2 Wg 



= 1 +T 



¥(£, o ,0)+TA g ¥(£„0) 
1 



+ 0(T 
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2^ 
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AV 
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+ 0(T 



AV 



(1 ^ 6) l+TAl + 0(T 2 ) 



(1.136) 
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where |3 is computed for A = 2. This calculation is based on the assump- 
tion that the expansion of the exponential commutes with the functional 
integration over geometries. Then we have: 

dim[T] = dim[V 2/dH ] = dim[V- p/<x ], (1.137) 

and thus 



d H = -^=2^Z^ + ^Z?, (1.138) 
(3 V25 - D + v/T^D V ' 

for matter with central charge D coupled to two-dimensional gravity. For 
D = 0, that means in the case of pure gravity, (1.138) predicts dn = 4, in 
agreement with (1.41). In the semiclassical limit D — > -co we expect dn = 2, 
which is also correctly predicted by (1.138). 

The validity of (1.138) has been widely discussed. The main assumption 
of the derivation sketched above is that one is allowed to interchange the 
expansion of e JAg with the functional integration in (1.136). Further input 
is the scaling assumption for (R 2 (T)) y . In principle the same result for the 
intrinsic Hausdorff dimension should be obtained if one uses some power of 
the Laplace operator. However this is not true and results obtained from 
higher orders of the expansion of e TAg contradict the result from the first 
order. 

An alternative theoretical prediction based on the transfer matrix method 
has been given in [34, 48, 49]. They get: 

2 24 

d H = — = =• (1.139) 

7 1 — D + a/(25 — D)(l — D) V ; 

For D = this yields again d H = 4, while for D — > — oo (1.139) predicts 
d H = 0. 

To settle this question, the authors of [50, 51] performed a high precision 
Monte Carlo analysis for D = —2. Their result is d H = 3.574(8), in agreement 
with the value dn = 3.562 from (1.138) and in clear disagreement to dn = 2 
from (1.139). This clearly rules out the validity of (1.139) for matter with 
D < coupled to gravity. 

However, for unitary matter with D > coupled to gravity the situation 
is less clear. Most simulations report values d H ~ 4 [52, 39, 53], and a back- 
reaction of the matter on the structure of the quantum spacetime cannot be 
observed. However, values extracted from the scaling of the matter correla- 
tion functions are higher, though less reliable. All data are consistent with 
the conjecture [52, 39]: 

d H = 4, for < D < 1, (1.140) 
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D = -5 


D = -2 


D =0 


° = i 


D = j 


D = 1 


ref. 


3.236 


3.562 


4 


4.212 


4.421 


4.828 


(1.138) 


1.236 


2 


4 


6 


10 


oo 


(1.139) 


3.36(4) 


3.574(8) 


4.05(15) 


4.11(10) 


4.01(9) 


3.8-4.0 


[54, 55] 








3.96-4.38 


3.97-4.39 




[55] 



Table 1.1: Current status of numerical values for the intrinsic Hausdorff dimension 
dn compared to the two different theoretical predictions (1.138) and (1.139). The 
values in the first part of the table are the theoretical predictions. In the second 
part of the table we summarize results from the scaling of the two-point function for 
the random surfaces, while in the last part the Hausdorff dimension is determined 
from the scaling of the matter correlation functions. See the references for technical 
details about the simulations and the data analysis. 

although the validity of (1.138) for D > is not completely ruled out. The 
current status of the numerical simulations is summarized in table 1.1. 



Chapter 2 

Discretization of 2d quantum 
gravity 

The formulation of quantum field theory in terms of renormalized Euclidean 
functional integrals leads to an identification of quantum field theory with 
statistical mechanics. A key ingredient for this identification is the discretiza- 
tion of spacetime. Then the powerful techniques from the theory of critical 
phenomena get available and have proved to be invaluable for the under- 
standing of renormalization and of non-perturbative phenomena. 

Therefore it is natural to attempt a discretization of the continuum the- 
ory of quantum gravity. However, spacetime no longer plays the role of a 
mere background but is the dynamical variable itself. This, combined with 
the diffeomorphism invariance of the continuum theory poses a problem for 
the discretization which has successfully been solved by the method of dy- 
namical triangulations [56, 57, 58]. The main idea is to discretize geometry 
directly, with no reference to coordinate parametrizations. That means that 
the functional integral over equivalence classes of metrics on a manifold M is 
replaced by a finite sum over piecewise linear spaces which are constructed 
by successively gluing d-simplices together. The fixed edge length of the 
simplices introduces a reparametrization invariant cutoff into the theory. In 
this way quantum gravity can be formulated as an ensemble of discrete ran- 
dom manifolds and the machinery of critical phenomena can be applied. The 
continuum theory is recovered at a critical point of this ensemble. Masses 
and continuum coupling constants are defined by the approach to the critical 
point in the scaling limit. 

Research in the formulation of classical gravity in purely geometrical 
terms was initiated in 1961 by T. Regge [59], while the first attempt at 
a discretization of quantum gravity in the sense of this chapter has been 



32 



2.1. DYNAMICAL TR1 ANGULATION 



33 



done by Weingarten in 1980 [60, 61]. A good review with some comments 
about the historical development of the subject is [62]. Other reviews can be 
found in [63, 64, 65]. Much of the material presented in this chapter is taken 
from these articles. In [66] the relation between the matrix-model technique 
and dynamical triangulation is reviewed. 

In the first section of this chapter we introduce the method of dynami- 
cal triangulations and discuss the scaling properties of the theory. Branched 
polymers, which have been mentioned before, are discussed in some detail. In 
section 2.2 we give a brief review of the matrix-model techniques which have 
been applied to solve two-dimensional quantum gravity. For more details 
and further issues we refer to the review articles given above. The explicit 
solution of discretized two-dimensional quantum gravity shows, that the con- 
tinuum theory is the scaling limit of dynamical triangulation. Therefore both 
theories can be identified. 

In section 2.3 we outline how the fractal structure and the scaling prop- 
erties of pure two-dimensional quantum gravity can be obtained. 

An alternative method for the discretization of quantum gravity, known as 
quantum Regge calculus, has been suggested. This is discussed in chapter 3 
of this thesis. 

2.1 Dynamical triangulation 

Dynamical triangulation is a discretization of quantum geometries. No ref- 
erence to parametrizations has to be made. At the same time it provides 
a regularization of the continuum theory by introducing an explicit cutoff. 
In this section we can only outline the beauty and power of dynamical tri- 
angulations in two dimensions. We begin with a short introduction into 
discretization of continuum manifolds and geometries. Then we define the 
partition function of two-dimensional dynamical triangulation and discuss its 
scaling properties. The susceptibility exponent y is characterized through 
the fractal structure of the random spacetimes. We end this section with 
a discussion of the scaling properties of the two-point function which are 
illustrated by means of the model of branched polymers. 

2.1.1 Discretization of geometry- 
Let us define an oriented n-simplex as an (n + 1 )-tuple of points modulo 
even permutations. Thus a 0-simplex is a point, a 1 -simplex is a pair of 
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points which can be identified with a line, a 2-simplex is a triangle, a 3- 
simplex a tetrahedron etc. A subsimplex a' of a simplex a is defined as 
a subset of cr: cr' C o\ By gluing simplices together one gets a simplicial 
complex. Simplicial complexes have no fixed dimension, that means they 
are not manifolds. However, their structure is sufficient to define exterior 
calculus on them, which physically means that scalar fields, gauge fields, 
antisymmetric tensor fields etc can be defined on them. 

A d-dimensional simplicial manifold is a simplicial complex such that 
the neighbourhood of each point in the complex is homeomorphic to a d- 
dimensional ball. In one and two dimensions each simplicial manifold can 
be obtained by gluing pairs of d-simplices along some of their (d — 1 )-edges. 
In higher dimensions this successive gluing yields in some cases only pseudo- 
manifolds, in which a neighbourhood of a point can also be homeomorphic 
to a topologically more complicated object. 

The structure of simplicial manifolds corresponds to a discretization of 
the continuous manifold structure. To define geometrical concepts we have 
to introduce a metric. A d-dimensional simplicial manifold can be equipped 
canonically with a Riemannian metric by demanding that the metric is flat 
inside each d-simplex, continuous when the (d— 1 )-subsimplices are crossed 
and that each (d— 1 )-subsimplex is a linear flat subspace of this simplex. 
A simplicial manifold equipped with such a metric is called piecewise linear 
space. 

The canonical metric can be defined by assigning lengths to the Ni 
links (1-subsimplices) of the simplicial manifold. First note that any d- 
dimensional simplicial manifold M can be covered with charts (U, 40 such 
that each d-simplex is parametrized by barycentric coordinates. Here 

U=RGR^ + - + ^<1}, (2.1) 

and 4> : U — > M is given by 

=£, 1 Vi + ... + ^d+(1 -^-...-WVd+i, (2-2) 

where y-\, . . . ,t)d+i are the coordinates of the vertices of the simplicial com- 
plex which live in some ambient space R n . Then on each chart the canonical 
metric is defined by 

This metric is Euclidean inside each d-simplex and continuous if a (d — 1 )- 
face is crossed. It is compatible with the manifold structure and can be 
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expressed solely in terms of the link lengths li, . . . , In v For a triangle (2.3) 
can be written as 

The intrinsic curvature of a two-dimensional piecewise linear space is 
located at the vertices, that means its O-subsimplices. The curvature can 
clearly not reside in the triangles, since they are flat. Furthermore the met- 
ric is continuous when crossing the edges of the triangles. Since we define 
curvature as an intrinsic concept of the piecewise linear space it should be 
independent of the embedding. We can bend the surface around an edge 
without changing the geometry - an observation which the German math- 
ematician C. F. Gauss made in his famous paper Disquisitiones generates 
circa superficies curva (1828) and which he called theorema egregrium. 
Therefore it is intuitive that the curvature cannot reside in the edges. Rather, 
it is concentrated at the vertices or in general at the (d — 2)-hinges of the 
piecewise linear space. To each vertex v in the surface we can assign a deficit 
angle 5 V as the difference between 2n and the sum of the angles meeting at 
v: 

6 V = In - a v>t , (2.5) 

trvet 

where the sum goes over all triangles t containing v. <Xv )t is the angle at v 
inside t. The scalar curvature TZ V at v is defined as 

v tvet 

where A t is the area of the triangle t. The scalar curvature 1Z in two dimen- 
sions equals two times the Gaussian curvature. With these definitions we 
finally have: 

N v 

total area of the surface = ^ dA v , (2.7) 

v=1 

and the Gauss-Bonnet theorem is now valid in the form 

N v 

Y_ ^ v dA v = 4ttx, (2.8) 

v=1 

which is equivalent to the polyeder formula N t — Nt + N v = x by Descartes 
and Euler. 
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2.1.2 Definition of the model 

The original proposal of Regge [59] consists in constructing a sequence of 
piecewise linear spaces to approximate a given smooth manifold M such 
that for any sufficiently smooth function f 



dA(£J f(£,), (2.9) 

M 



where f v is the value of the function f at the vertex v. The applicability of 
this point of view in the quantum case will be discussed in chapter 3. Here 
we take a somewhat different point of view. Our aim is to integrate over 
all equivalence classes of metrics on a given manifold M. We introduce a 
reparametrization invariant cutoff by assigning the length a to each edge of 
piecewise linear spaces. We construct all possible piecewise linear spaces by 
gluing equilateral triangles together. In this case the formulas above simplify 
considerably. Let n v be the order of the vertex v. Then we have: 

\H> In 
dA v = -^a 2 n v , ^ v dA v = — (6 - nv). (2.10) 

For simplicity we rescale a 2 by ^ and then set a = 1 . (2.10) shows that two 
triangulations which cannot be mapped onto each other by a relabelling of 
the vertices lead to different local curvature assignments. Thus they define 
different metric structures. The set of combinatorially non-equivalent tri- 
angulations defines a grid in the space of diffeomorphism classes of metrics. 
The proposal described in this chapter relies on the hope that this grid gets 
uniformly dense in the limit a — > 0. 

To define a regularized theory of quantum gravity we replace the action 
S(g,A) = AV g in (1.6) by 

S T (n) = jiN t , (2.11) 

where |x is the bare cosmological constant. The integration over equivalence 
classes of metrics on M is replaced by a summation over all non-isomorphic 
equilateral triangulations T with the topology of M. In the case of mat- 
ter coupled to quantum gravity, the matter fields X can be defined on the 
vertices, the links or on the triangles. In general the action S matte r and the 
functional integration over the matter fields will depend on the triangulation 
T. The discretized version of the partition function can thus be written as 

Z{\i) = Y -Le~ STM \v T [X] e - s -«" (T - x) , (2.12) 
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where Cj denotes the symmetry factor of the triangulation T which is equal 
to the order of the automorphism group of T. The summation goes over a 
suitable class T of abstract triangulations T defined by their vertices and a 
connectivity matrix. Appealing to universality details of the chosen class of 
triangulations such as whether closed two-loops are allowed or not should 
not be important for the theory. This has been verified a posteriori. 

For simplicity we will ignore any possible matter fields coupled to gravity 
in the following equations. 

Then we can write: 

oo 1 

Z(^=2>-^ N Z(N), Z(N)=£— , (2.13) 

N=l TeT N T 

where 7^ denotes the subset of triangulations with N triangles of T. Expec- 
tation values of observables O are defined as: 

(O)^ = ^^L^-^Od), (2.14) 

( >N = zk)Z- h 0[Tl (2 - 15) 

TeT N 

in the "grand-canonical" ensemble with fixed cosmological constant and in 
the "canonical" ensemble with fixed volume respectively. \i can be under- 
stood as the chemical potential for adding new triangles to the surfaces. 

Z(N) can be interpreted as the number of triangulations in T N . This 
number is exponentially bounded: 

Z(N) = e^ N N y - 3 (1 +0(1/N)). (2.16) 

The proof for spherical topology has been given in [67] . For a general proof 
by purely combinatorial methods in the spirit of Tutte see [62]. It can be 
proved that the critical point \x c does not depend on the topology. 

That means that the statistical ensemble defined by (2.12) has a critical 
point [L c . Z(|x) is analytical for \jl > \i c and contains non-analytical parts 
at the critical point. The latter have our biggest interest since they are the 
universal parts given by the large N behaviour of Z(N). The continuum limit 
should be defined as \i — > \x c . Close to the critical point we have: 

Z((x) = [i c ) 2 ~ y + less singular terms, ( 2 -l?) 



from a discrete Laplace transformation of (2.16). Many analytical and nu- 
merical investigations have revealed that y assumes the same values as in the 
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continuum theory of two-dimensional quantum gravity and is given by (1.14). 
Thus the susceptibility exponent y can also be characterized as the sublead- 
ing power correction of the exponentially growing number of triangulations 
for a fixed topology. Furthermore, y can be characterized by the branching 
ratio into minimal bottleneck baby universes [68] , that means parts of the 
universe which are connected to the rest by a three-loop of links which do 
not form a triangle, see figure 2.1. The smaller part is then called minimal 




Figure 2.1: Branching of a minimal bottleneck baby universe from the 
larger parent. 

bottleneck baby universe and the larger part is called parent. In the ensemble 
of universes with volume N the average number (JV(V)) N of baby universes 
with volume V is given by 

(.W0) N ~~ ^qyVZ(V) (N - V)Z(N - V). (2.18) 

The factors Z( V) and Z(N — V) are the weights of the baby universe and of 
the parent respectively. The volume factors V and N — V reflect the fact that 
the baby universe can be attached at any triangle. Finally, 3! is the number 
of ways the two boundaries of the parent and of the minimal bottleneck baby 
universe can be glued together. (2.18) is valid in the generic case where there 
are no additional symmetries. Assuming that Z(N) is given by (2.16) we get: 

(W0> N - NV^- 2 (l - XV 2 . (2.19) 

For V « N this reduces to (A/"(V)) N ~ NV Y ~ 2 . It is interesting to note 
that y is a function of the matter coupled to gravity. Thus it describes an 
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aspect of the fractal structure of quantum spacetime that originates from the 
back-reaction of the matter on gravity. 

2.1.3 The two-point function 

In analogy to the continuum formalism in chapter 1 we can define a two- 
point function which in a natural way contains information about the fractal 
structure of the random surfaces. 

In principle, we could define a geodesic distance by using the canonical 
metric described in section 2.1.1. Alternatively, the use of simplified defini- 
tions has been suggested. The geodesic distance between two links li and I2 
is denned as the length of the shortest path of triangles connecting li and 
I2. The geodesic distance between a link li and a set of links C is defined 
as the minimum of the geodesic distances between li and the elements of C. 
Furthermore we define the geodesic distance between a loop C\ and a loop 
£2 to be r if all links of C-\ lie a geodesic distance r from the loop C 2 . Note 
that this definition is not symmetric in C\ and C 2 . Similarly we could have 
defined the geodesic distance as the distance between vertices or between 
triangles. 

Then the two-point function is defined as: 

00 

G(mr) = ^ e^ Nt = ^e^ N G(N,r), (2.20) 

TeT(2,r) N=l 

G(N,r) = 1, (2.21) 

TST N (2,r) 

where T(2,r) denotes the class of triangulations T with two marked links 
separated a geodesic distance r, and 7n(2,t) the subclass where all triangu- 
lations consist of N triangles. Note that r is an integer in units of the lattice 
spacing. G{\i,r) is the discretized form of G(A,R) in (1.23), while G(N,r) is 
the discretized form of G(V, R) in (1.25). 

An important property of the two-point function is: 

G[[i,r-t + r 2 ) > const x G[[i,r-t) G[[i,r^, (2.22) 

which follows from the fact that each term on the right hand side uniquely 
defines a term on the left hand side. This is demonstrated in figure 2.2. From 
this one can conclude that the limit 



-logG(M.,r) 
hm = m((j.J 

T — > OO X 



(2.23) 
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Figure 2.2: Two surfaces with two marked links separated a geodesic 
distance ri and rz respectively can be glued together after cutting open 
one of the marked links each to form a surface with two marked links 
separated a geodesic distance ri + rj. 



exists with m(u-) > and m/(u.) > 0, for \i > \i c . Thus the two-point 
function falls off exponentially for r — > oo. These relations are by now well 
known. However, the complete proofs require some technical arguments, 
compare [62]. 

Below we will show that a continuum limit of the discretized theory can 
only exist if m(|j.) scales to zero as \i — > u c . Let us assume that this is the 
case: 

m(u) - {[i-\i c y, for \i c . (2.24) 

In general we expect that close to the critical point the exponential decay of 
G(fx, r) turns over into a power law. More precisely: 

G(u,r) - e - mMr , for 1 < m(u)r, (2.25) 
G(mr) - r 1_1 \ for m(ji) < m{\i)r < 1. (2.26) 

The exponent r| is called the anomalous dimension in analogy to ordinary 
statistical systems. Furthermore we can define the susceptibility x(M-) in the 
discretized ensemble as 

oo 

xM = X G ^> r )> ( 2 - 27 ) 

r=1 

Close to the critical point, where triangulations with large N t dominate and 
symmetry factors play no role, we have 

Xlv)-^=^-(v-\Lc)- y . (2.28) 
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Similar to chapter 1 we can introduce several concepts of fractal dimensions 
to characterize the geometrical structure of the two-dimensional quantum 
spacetime. The Hausdorff dimension dn in the grand-canonical ensemble is 
defined by 

(N) T ~ r dH , for r -> oo and m{\i)r = const, (2.29) 
where the average (N) T is given by 

M^rT-i L N t e-^=- ai °« G( ^ r) . (2.30) 
Gfu.rj ^— Oil 

^' ' TeT(2,r) ^ 

If we perform the derivative under the constraint in (2.29) it follows that 

(N) r ~mV)r~r^. (2.31) 
Thus the Hausdorff dimension d H is related to the scaling exponent -v by: 

y = ^-. (2-32) 
With (2.28) we can derive Fisher's scaling relation 

Y = v(2-n), (2.33) 

which relates the critical exponents denned above. From the long distance 
behaviour of G((x,r) we can compute the long distance behaviour of G(N,r) 
by a saddle point calculation: 

G(N,r) -e^^^, (2.34) 

with c = — 1 ) "v"^ . 

Analogously to the continuum formalism in chapter 1 we can define a 
Hausdorff dimension dn by the short distance scaling of the two-point func- 
tion in the canonical ensemble with fixed number N of triangles. For r = 0, 
G(N,0) is the one-point function and behaves as 

G(N,0) ~ e McN N Y ~ 2 . (2.35) 

This is because for large N the one-point function is proportional to NZ(N) 
since it counts triangulations with one marked link. 

Now let r <^ N 1/dk and count the number n(r) of triangles which lie a 
geodesic distance r away from a marked link. The average of n(r) in the 
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ensemble of surfaces with one marked link defines the Hausdorff dimension 
d h by 

(n(r)) N ~^|^~A- 1 ) for1«r«N 1 ^. (2 . 36) 

The first relation follows from the definition of (n(r)) N . Together with (2.35) 
we conclude: 

G(N,r) ~ r dH_1 N T_2 e vlcN , for1<r<N 1/d \ (2.37) 

Finally, the short distance behaviour of G((x,r) follows by a discrete Laplace 
transformation. Close to the critical point we get: 

G(mr) -r^- 1 . (2.38) 

Together with (2.26) this proves again Fisher's scaling relation. Furthermore 
it turns out that both definitions of the intrinsic Hausdorff dimension are 
equivalent in two-dimensional quantum gravity. However, counter examples 
exist, as will be demonstrated in the next section. 

The mass m( \x) determines the scaling of the discretized theory. To see 
how continuum expressions can be approached in the discretized theory let 
us reintroduce dimensions. The renormalized cosmological constant A in 
terms of the bare cosmological constant \i is defined by: 

(j. — \x c = Aa 2 , that means a{\x) ~ (J. c ) 1/2 - (2.39) 

If the mass shall survive in the continuum limit, it has to be introduced as: 

M = m(n) a{[L)- 2 ^ = cA v , (2.40) 

with a constant c. If the two-point function shall survive, we must have 

m(|i)T = MR, (2.41) 

where the continuum parameter M and R are held fixed and the number of 
steps r goes to infinity, that means: 

R = Ta([ip. (2.42) 

From (2.25) and (2.26) we see that the continuum two-point function is given 
by: 

G(A, R) = lim a((j.) 2 ^ (1_Tl) G((j.,T) ) m(n)r = MR. (2.43) 

This scaling form can be explicitly verified by an exact calculation of G((x,r) 
for pure gravity [22, 23], compare section 2.3. 
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2.1.4 Branched polymers 

An interesting simple example for the scenario described above is provided 
by the model of branched polymers which is important in the understanding 
of the phases of two-dimensional quantum gravity. The critical indices have 
been computed in detail in [69, 21]. The partition function is defined as 

Zt^Ipe^ 1 n f ^> ( 2 - 44 ) 

BP ar veBP 

where the sum goes over all branched polymers, that means all connected 
planar tree graphs. Nv denotes the number of links in a branched polymer. 
The product goes over all vertices in the tree graphs and n v is the number 
of links joining at the vertex v, called the order of v. Usually, the (unnor- 
malized) branching weight f (n) is a non-negative function of the order of 
the vertices. Finally the symmetry factor Cbp is chosen such that rooted 
branched polymers, that means with the first link marked, are counted only 
once. 

The one-point function G{\i) is defined similarly as the sum over the 
rooted branched polymers. In this case, the symmetry factor drops out. 




Figure 2.3: Illustration of the self-consistent equation for rooted 
branched polymers. The first link gives a factor e~^. At the vertices 
the branching weight has to be included. 



From figure 2.3 we follow that G(|x) satisfies the equation 

G(u] = e-*(1 + f(2)G(n) +f(3)G(^) 2 + ...) 
We can solve this relation for as a function of G((x): 
^ 1 +f(2)G(^+f(3)G(^) 2 + ... _ f(G) _ 



= HG). 



(2.45) 



(2.46) 
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G(Hc) G(n) 



Figure 2.4: Graphical solution of (2.46). The critical point is identified 
as the minimum of as a function of G. 

This equation is illustrated in figure 2.4. The critical point is at the minimum 
of F(G). Since all f(n) are positive the minimum is unique and satisfies 
F'(G(Hc)) = and F"(G(n c )) > 0. Therefore we have: 

G{\i)~G[\i c )-c{\i-vl c )?, for^^c, (2.47) 

where c is some constant. Since G((x) is the one-point function we would 
have expected a behaviour G[\i) ~ M. C ) 1_Y - Thus we conclude that the 
generic value of y for branched polymers is \. 

If we allow some of the values f (n) to be negative, we can fine tune the 
f(n) such that the minimum of F(G) satisfies 

F (lc) (G(^ c )) =0, for k= 1,... ,m- 1, and F (m) (G(^)) ^ 0. (2.48) 

This model is called m'th multicritical branched polymer model. For m = 
2 we simply recover the ordinary branched polymer model. The critical 
behaviour of the one-point function is changed to 

G[\i) - G(mO -cfji-jic)™, for ^^ jig, (2.49) 

and we get y = 

To define the two-point function G(fx, r) for branched polymers, let the 
geodesic distance between two marked vertices x and y be the minimal length 
of paths between the two points. This is of course unique. The two-point 
function is defined as the partition function of branched polymers with the 
constraint that there are two marked points separated a geodesic distance r. 
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The path between the marked points can be viewed as a random walk where 
at each vertex a rooted branched polymer can grow out. If we think only in 
terms of intrinsic geometrical properties this leads to [21]: 

G(mr) ~ (e-^'(G(^)) r = e - mM \ (2.50) 
from which we conclude, using e _ ^jF'(G((x)) = 1 — that 

Tn((x) = k(h- \Lc)^, (2.51) 

in the general case for m > 2, with a positive constant k. Therefore it follows 
with (2.32) that the intrinsic Hausdorff dimension dn equals -^j. For m = 2 
we have d H = 2 and the branched polymers have the same intrinsic dimension 
as smooth two-dimensional manifolds. For m — > oo, dn approaches one, in 
agreement with the fact that in this limit the branched polymers approach 
ordinary random walks. The two-point function G[[i,r) is given by 

G(mr) = const x e - Kr{ ^ ' ^ . (2.52) 

We can also compute the canonical intrinsic Hausdorff dimension d h for 
branched polymers defined in the ensemble of branched polymers with a fixed 
number N of links. The volume (n(r)) N of a spherical shell of geodesic radius 
r is defined in analogy to (1.35) and (2.36). Z(N) for branched polymers is 
given by an inverse Laplace transformation of (2.44) and scales as Z(N) ~ 
Similarly the two-point function G(N,r) for graphs with fixed 
volume N is given by an inverse Laplace transformation of G((x, r). Inserting 
the scaling behaviour and taking only the leading orders we find for small 
r [51]: 

(™M)n r-r. for N ^ < r < N 1 "^- ( 2 - 53 ) 

This shows that d^ = 2 for all values of m. Thus for the generic branched 
polymers (m = 2) the canonical and the grand-canonical definition of the 
intrinsic Hausdorff dimension give the same result. However, this is not 
true for the m'th multicritical branched polymer model with m > 2, where 
d H = and d H = 2. 

The fractal structure of branched polymers embedded in M° is described 
by the extrinsic Hausdorff dimension D H . D H can be computed from the 
scaling of the embedded two-point function. The path between two marked 
points x and y is an embedded random walk at which at each point branched 
polymers can grow out. The integration over the embedding field is the same 
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as for the ordinary random walk in D dimensions. The branching at the 
vertices corresponds to a renormalization of e~^. Prom this renormalization 
D H can be computed with the result D H = | = ^ [21]. Thus D H = 4 for 
ordinary branched polymers and D H — > 2 for m — > oo. 



Major advances in the theory of two-dimensional quantum gravity have been 
made when it was realized, that the theory is integrable in the discrete for- 
mulation, even for some types of matter coupled to the surfaces [56, 70, 57, 
71, 72, 58, 73]. 

In the development of these results the formulation of dynamical trian- 
gulation in terms of matrix models has been of some importance. The Her- 
mitean one-matrix model has a direct interpretation in terms of randomly 
triangulated surfaces. This is demonstrated in the beginning of this section. 
Also the Hartle-Hawking wavefunctionals have a direct natural formulation 
via matrix models. The Dyson- Schwinger equations for these models allow 
the treatment of the theory by methods of complex analysis. These so-called 
loop equations can alternatively be derived using only combinatorial meth- 
ods [62] . The explicit solution of the loop equations provides deep insight 
into the properties of the theory and allows the computation of the scaling 
limit of the Hartle-Hawking wavefunctionals. 

We cannot review the vast literature on numerical simulations in quantum 
gravity in this work. However, let us remark that all numerical simulations 
of two-dimensional dynamical triangulation are consistent with the scaling 
hypotheses, with universality, and - most remarkably - with the continuum 
theory of two-dimensional quantum gravity. 

2.2.1 Dynamical triangulation by matrix models 

The discretized partition function of two-dimensional quantum gravity can 
be represented as an integral over a Hermitean N x N -matrix Consider 
the Gaussian integral 



2.2 Matrix models 




(2.54) 



where the measure deb is denned as 



d4)=n d ^n d:5c ^- 



(2.55) 



i<3 i<j 
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The propagator of (J) is given by 

<<h*tT> - 



d$> e^* 2 <bij(b V y = Sij/Sjv, (2.56) 



with Z(N , 0) = Jdcf) e _ ? Tr * 2 . Diagrammatically, the propagator can be repre- 
sented as a double line where the two lines are oriented in opposite directions. 
The integral (2.54) is performed by doing all possible contractions of the K 
Trcf) 3 vertices. The dual picture of this corresponds to gluing K triangles 
together to form all possible closed and not necessarily connected surfaces 
of arbitrary topology. The contribution from a particular graph forming a 
closed surface is N Nv , since each vertex contributes with a factor N. If we 
make the substitution 

Tr4> 3 -> -^=Trct) 3 , (2.57) 



each closed surface with Euler characteristic x contributes with g K i\j N v-K/2 = 
g K N x . We get the sum over arbitrary closed surfaces with any number of 
triangles if we sum (2.54) over K: 



z(N, g ; 



dcM-^* 2 ^^* 3 . (2.58) 



Taking the logarithm of this projects on the connected surfaces only. Thus 
we see that 

with 

1 = log N , and \i = - log g (2.60) 

is a formal definition of the discretized partition function of two-dimensional 
quantum gravity, including the sum over topologies. Here G is the gravita- 
tional coupling constant. The matrix integral (2.58) is of course not conver- 
gent, but it has been suggested, that a closed form like (2.58) might define 
a non-perturbative definition of the sum over topologies after analytic con- 
tinuation. In general however, the integral will be complex [74, 75] and the 
problem of summing over topologies has not yet been solved. 

Equation (2.58) admits a l/N 2 -expansion which is identical to an ex- 
pansion over topologies. Here we will restrict ourselves to spherical topology 
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which corresponds to the leading order in 1 /N 2 . It is convenient to generalize 
the matrix integral (2.58) to 



Z(N, 9l ) = 



dcp e- 1 ^*), (2.61) 



V(4>) = X— * n > ( 2 - 62 ) 

Tl=1 

where we have performed the rescaling cf) — > a/N^- (2.61) is interpreted as 
a perturbative expansion around a Gaussian integral. Thus the Gaussian 
coupling constant is bigger than zero, while the others are chosen smaller 
than zero. The expectation value of some observable 0{§) is denned as 

1 



(0(<M> 



d* e^^O^). (2.63) 



Z(N, 9i ) 

Differentiating the logarithm of (2.61) with respect to the coupling con- 
stants g n defines expectation values of observables like Tr(b kl • • • Trct) kb . They 
are interpreted as the sum over all surfaces which have b polygons with lq 
links as their boundary. The generating function of the connected expecta- 
tion values of these observables is given by 

k,,...,lc b =0 Z 1 "' Z b 

= N b - 2 /Tr— '-T---TT — ] —) . (2.64) 



' z, - 4> z b - 4> , 

In the large-N limit this is the generating function for the discretized Hartle- 
Hawking wavefunctionals, whose continuum counterparts have been defined 
in section 1.1.2. In fact, already the one-loop correlator W(z) contains all 
necessary information, since the higher correlators can be computed from it 
by differentiation: 

w(z ' ^ = dvL)dV(L '--d^) W(Z|) ' (2 ' 65) 

where the loop insertion operator is defined as 



dV(z) ~~ ^ z n+1 3g 



Let us define the density p(A) of eigenvalues of the matrix integral (2.61) by 



N 



p(A) = (^6(A i -A)), (2.67) 



i=1 
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where the At are the N eigenvalues of the matrix cp. Then the one-loop 
correlator W(z) can be written as: 



W(z) 



dA (2.68) 



For polynomial potentials the support of p(A) is confined to one or several 
intervals on the real axes in the limit N — > oo. The solution in the case of 
several cuts has been given recently [76] . Here we will only deal with a single 
cut [a, b] on the real axes. Then W(z) will be analytic in the complex plane 
with a cut at the support [a, b] of p(A). It follows that 

2m p(A) = lim (W( A + ie) - W(A - ie) ) . (2.69) 
2.2.2 The loop equations 

The matrix model (2.61) can be solved by many methods. It is most conve- 
niently and systematically done by using loop equations, the Dyson-Schwinger 
equations for the matrix models. Let us consider the transformation 

cp^cp + e— !— ■ , (2.70) 

z — cp 

which makes sense if z is not an eigenvalue of cp and real, so that the new 
matrix remains Hermitean. Then the measure and the action transform as 

dcj) -> (1 + eTr— !— Tr— *— ) dcp, (2.71) 

V z-4> z-$J 

TrV(cp) -> TrV(cp) + eTr^-^, (2.72) 

z — cp 

to first order in e. Since the matrix integral (2.61) is invariant under such a 
definition of the integration variable we get: 

d(J) e -mw ( *) |^ Tr _J_y_ NTr YMj =0- (2.73) 

The second term can be rewritten as a contour integral involving the one- 
loop correlator by using the eigenvalue density p(A), while the first term is 
related to the two-loop correlator. This leads to the standard form of the 
loop equations [77]: 

I^!M W ( t u)=W( Z ) 2 4w(z,z). (2.74) 
rlnx z — iv N z 



50 



CHAPTER 2. DISCRETIZATION OF 2D QUANTUM GRAVITY 



The contour C goes around the cut [a, b] but does not enclose z. The loop 
equation can be solved systematically to all orders in 1/N 2 , that means for 
all genera of the random surfaces [78, 79]. For spherical topology, that means 
in the limit N — > oo, (2.74) simplifies to 

o^XM Wo (cu)=W (^, (2.75) 
t c 2n\ z — tu 

where the subscript denotes the genus of the surfaces. By the definition 
we know that W (z) = \ + 0(z~ 2 ). Therefore (2.75) can be solved by a 
deformation of the contour C to infinity. For a polynomial action V of degree 
n we get: 

W (z) = 1 (v'(z) - V(z - a) (z - b) f_ M k [a, b, gj (z - b)^ 1 j 

f duo V'{w) J{z-a){z-b) 
= o — 1 J y_ (2.76) 
J c 2m z-a) ^(co-ojfco-b) 

where the moments M k are defined as 

M k [a, b, gj = i ^ V ' [a)) (2.77) 

J c 2m (cu _ a )i (cu - b) v+ 2 

and vanish for k > n— 1 . The endpoints a and b of the cut are self-consistently 
determined by the equations 

M_i[a,b,gd = 2, (2.78) 
M [a,b, gi ] = 0, (2.79) 

which are a consequence of (2.76) together with W (z) = \ + 0(z~ 2 ). 

For the matrix model (2.58) which directly corresponds to triangulated 
surfaces, the solution for W (z) is: 

W (z) = lv'(z)+f(mz), 



f(mz) = |(z-c)V(z-a)(z-b), (2.80) 
where a = e _,x and c = - — 

» g 2 

In principle one could compute all higher loop correlators from (2.76) by 
applying the loop insertion operator (2.66), which gives a complete solution 
of two-dimensional dynamical triangulation for spherical topology. 
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In fact, a closed expression for the b-loop correlators can be obtained if 
the small scale details of the theory are adjusted. Instead of the Hermitean 
matrix model one considers the complex matrix model with a potential 



n=1 



TrV(cb+ci)) = V ^Tr(ct) + cb) n , (2.81) 
and an integration measure 



N 



deb = ] J dD^cbijdJcby. (2.82) 

i,j=1 



The one-loop correlator is now defined as 



n=0 

higher loop correlators are denned analogously. The term Tr(cb + 4)) n can 
be interpreted as a 2n-sided polygon whose links are alternately coloured 
black and white. These polygons can be glued together as in the Hermitean 
matrix model with the additional constraint that black links have to be glued 
to white links [80]. Such short distance details about the gluing should be 
unimportant in the continuum limit. 

The loop equations for this model have been derived in [81]. Because of 
the symmetry (J) — > — cj) we have a = — b for the cut of the one-loop correlator. 
The solution for spherical topology is given by [81]: 

1 



W (z) = -(v'(z)-M(zWz 2 -b 2 ) t 



M z = o — — \ ' =YM k [b, gi \(z 2 -b 2 ^ \ 

J Coo 4m [cv 2 -z 2 )Vw 2 -b 2 

M k b,gJ = o— (2.84) 

J c 4m (co 2 -b 2 ) k+ ^ 

Here is a contour around the cut pushed to infinity. The position of the 
cut is determined by the equation 



M [b,gJ =2. 



(2.85) 
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The higher loop correlators are given by the following expressions [81]: 

w (zi,z 2 ) = T7 7 1 7 , 7 1 zum — h+ z um — h~ iz ^2 



W (Z 1 ,Z2,Z 3 ) 




16Mi y/(z$-b 2 )l4-b 2 ){zl-b 2 y 



(1 d \ ^ ^ 1 b 2 
Midb^j 2pM7n 2 (^- b 2)3/2 - ( 2 - 86 ) 

In these formulas all dependence on the coupling constants is hidden in Mi 
and b. Similar and only slightly more complicated statements are valid also 
for the Hermitean matrix model. 

2.2.3 Scaling limit and computation of y 

It has been discussed in section 2.1.3 that the ensemble of triangulated ran- 
dom surfaces has a critical point u c . The continuum limit with a renormal- 
ized cosmological constant A is approached by [i — > |~i c with u — u c = Aa 2 , 
where a is the lattice spacing. In terms of the coupling constant g for the 
model (2.58) this relation is g c — g ~ g c Aa 2 , close to the critical point. For 
the generalized matrix model (2.61) with n coupling constants we expect 
no qualitative changes except that the theory will be critical on a (n — 1 )- 
dimensional hypersurface. This hypersurface is identified by 

M^bteO.gJ =0, (2.87) 

since the b-loop correlators W {z-\, ... ,z b ) are exactly divergent for Mi = 0. 
Let us denote a point in the critical hypersurface by g C) i,i = 1 , . . . ,n and 
the corresponding endpoint of the eigenvalue distribution by b c . If we move 
slightly away from the critical surface 

9i = 9c,i(1 - Aa 2 ) = g Cii + 6gi, (2.88) 

there will be a corresponding change b 2 -> b 2 + 6b 2 which can be computed 
from (2.85) to be 

(Sb 2 ) 2 = -jr^ .Aa 2 ~ 6g t . (2.89) 

3M 2 [b c ,g C)i ] 

We rescale the cosmological constant such that b 2 = b 2 — av^. Because zi 
appears in (2.86) always in the combination (z 2 — b 2 ) it is natural to introduce 
a scaling of z t by 

z 2 = b 2 +aZ t . (2.90) 
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Since W (zi , . . . , z b ) is the generating function for b-loop correlators, we can 
compute the transition amplitude for b one-dimensional universes of lengths 
ui , . . . , n b from it by a multiple contour integration. The physical lengths of 
the loops will be Li = n t a. We see that in the continuum limit the number 
of links Tii on the boundaries has to go to infinity as a — > if the loops shall 
survive. By following this procedure one gets the generating functional for 
macroscopic b-loop amplitudes [81]: 

W {z u ... ,z b ) ~ a 5 -^W (A,L 1( ... ,L b ), 
W (A,Zi,...,Z b ) = (--^ -^=fl(Z^ + VX)-^ (2.91) 

for b > 3. The inverse Laplace transform of Wo (A, Z-\,. . . , Z b ) in the vari- 
ables Zi gives the Hartle-Hawking wavefunctionals of two-dimensional quan- 
tum gravity. It can be computed from (2.91) to be: 

W (A,L 1 ,...,L b ) = (-Ay 3 J_0^7[;e-^ + - +L ^. (2.92) 

Since the b-point function in two-dimensional quantum gravity should scale 
as A 2 ~ h ~ y , we can directly read off from this expression that 

y = -- (2.93) 

for pure gravity, in agreement with the KPZ-formula (1.14). These calcula- 
tions can be generalized to arbitrary topology [77, 78, 82, 79] and to models 
with matter coupled to gravity [83, 84, 62] by considering multicritical mod- 
els which are obtained by a fine-tuning of the critical coupling constants such 
that higher moments vanish as well as M|. 

It is a major result that all values of the susceptibility exponent y com- 
puted by these or other methods in the model of dynamically triangulated 
quantum gravity agree with the continuum formula (1.14) which has been 
derived in section 1.2.3. In fact, all calculations which can be done by dynam- 
ical triangulation and in the continuum approach have so far yielded the same 
results. Therefore we identify both theories as the theory of two-dimensional 
quantum gravity. Objects like the Hartle-Hawking wavefunctionals and other 
correlation functions can much easier be obtained in the discretized approach 
whose scaling limit yields the theory of two-dimensional quantum gravity. 
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2.3 The fractal structure of pure gravity 

Another major advance in two-dimensional quantum gravity was made by 
the explicit computation of the two-point function for pure quantum gravity 
by constructing a transfer matrix [22] and later by using a peeling decom- 
position [23]. The two-point function is a natural object since it provides 
all details about the scaling properties of two-dimensional quantum gravity. 
The astonishing result of these investigations is that the intrinsic Hausdorff 
dimension of two-dimensional quantum gravity is four. That means that 
even the dimensionality of spacetime is a dynamical quantity. 



2.3.1 The geodesic two-loop function 

Let us define with T(li , I2, r) the class of triangulations with an entrance loop 
li with one marked link and an exit loop I2, separated a geodesic distance r. 
We will also denote the number of links of li and l 2 with the same symbols. 
Then the geodesic two-loop function is denned as 

G(^,r;l 1 ,l 2 )= e-^\ (2.94) 

TeT(l, ,l 2) r) 

For r = we have the initial condition: 

G(M;l 1( l 2 ) =6i,,i 2 - (2-95) 

We introduce the generating function for G((x, r;l-|, l 2 ) by 

00 

G(^r;z 1 ,z 2 )= £ z^+V^G^r;^!,), (2.96) 
1, ,i 2 =i 

with the initial condition 

G(mO;zi,z 2 ) = - — f ( 2 - 97 ) 

Z]Z 2 Z]Z 2 — I 

By a two-fold contour integration the geodesic two-loop function can be 
reconstructed: 



dz-| 



c, lni 



dz 2 



G(^r;l 1 ,l 2 )=(|) ^4 2G (^; z i- z 2). (2-98) 



C z 



It is an important observation that the two-point function G ( fx, r) can be ob- 
tained from the geodesic two-loop function. Consider the two-loop function 
with an entry loop of length li = 1 which is equivalent to a marked link, 
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and an exit loop of arbitrary length I2 separated a geodesic distance r. We 
can close this surface by gluing a disc to the boundary \ 2 . The amplitude 
W (l2) of the disc can be computed from the one-loop correlator W (z) by 
a contour integration. An additional factor of I2 arises because we have to 
mark one of the I2 links of the exit loop. Thus we get: 



G(mr) = ^G(^r;1,l 2 ) l 2 W (l 2 ) 
i 2 =i 



= o 



dtu 1 

c 2m cu 



z 2 G(m.,t;z,— ) 
cu 



3cu 



cuWq(cu) 



(2.99) 



where C is a contour around zero. Thus we see that all information about the 
scaling of the theory can be obtained if the generating function of geodesic 
two-loop amplitudes can be computed. 

By a step- by-step decomposition of triangulations in T(li,l 2 ,r) a differ- 
ential equation for G((x, r;z, cu) can be obtained . It is intuitive that any tri- 
angulation in T(li , l 2 , r) can be decomposed into r rings of thickness one [22]. 
This leads to a transfer matrix formalism for the two-point function. Alter- 
natively, one can decompose the triangulations by a peeling procedure [23]. 
The differential equation one obtains is not exact but should be valid close 
to the critical point of the theory: 



3 9 

— G(mr,z,a>) = -2— [f(mz)G(mr,z,a>)] , 
where f(mz) is given by (2.80) for triangulated surfaces. 



(2.100) 



2.3.2 Scaling of the two-point function 

The differential equation (2.100) can be solved by the method of character- 
istic equations. The result is: 



G(m r; z,cu) 



f(n,2) zcu 



f(mz) z-cu-T 

where £ is the solution to the characteristic equation 

dz(z,r) 



dr 



given by 



1 



1 8, 



2f(m*), 



1 



£(z, r) c c sinh 2 ^_ §oT + Sinh -T ^/_Ji_ _ 52 j + 5. 



(2.101) 



(2.102) 



(2.103) 



56 



CHAPTER 2. DISCRETIZATION OF 2D QUANTUM GRAVITY 



G[\i,r) = const x 6o&i _ ; , 3 c (1 + O(6 )). (2.107) 



Here the positive constants 6i scale as 

6 = l y /( c - a ){c-b)=0{{\L-\i c )l) t (2.104) 

6, = (C "; )(C T b) =0(^-H,)i) t (2.105) 
c(b — a] 

62 = -^^ = 0(Ui-iic)i). (2.106) 

Now the two-point function G{[i,r) can be computed by inserting (2.101) 
into (2.99). Close to the critical point one obtains [85]: 

cosh(6 r) 
1 sinh 3 (6 r) 

Prom (2.104) we get: 

3 cosh (c[\i- |J.c)*r) 

G{yi,r) = const x (^- n. c )? ^ r ^- ) (2.108) 

sinh 3 f c(jx— M. c ) ? rJ 

where c = v^e - ^ is a nonuniversal constant. To read off the critical be- 
haviour of the theory we only have to analyze the asymptotic behaviour of 
the two-point function for large and for small distances. We see that: 

1 

• G{\i,r) falls off exponentially as e~ 2c( ^~^ )4T , for r — > 00. Thus the sus- 
ceptibility exponent equals -v = \ and the intrinsic Hausdorff dimension 
d H of pure two-dimensional quantum gravity equals four. 

• For 1 r (\i — |x c )~4 the two-point function behaves like r , 
that means the anomalous scaling dimension is r\ = 4. For ordinary 
statistical systems r| is smaller than two. It is remarkable, that the 
critical exponents of the theory still satisfy Fisher's scaling relation 

7 = y{2-r[). 

• Any geodesic two-loop function G((j., r; li, l 2 ), which can be computed 
from (2.98) and (2.101) fulfills the same scaling relations as the two- 
point function, provided that li and I2 stay finite in the limit \i — > \x. 



c- 



The continuum limit of the two-point function is given by 



3 cosh ( cA4 R 

G(A,R) = CA? ) 7 (2.109) 

sinh 3 ( CA4RJ 



where C is a constant. 
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2.4 Conclusion 

In this chapter we have demonstrated, how to use discretized systems to de- 
scribe two-dimensional quantum gravity. In two dimensions, the continuum 
theory of quantum gravity and the theory of dynamical triangulations agree 
with each other in that region, where the continuum theory can be evaluated, 
that means for matter with central charge D < 1 coupled to gravity. The 
continuum limit of dynamical triangulations in two dimensions can therefore 
be identified with the continuum theory. 

In fact, the discretized approach is more powerful than the continuum 
approach in the sense that natural observables like the Hartle-Hawking wave- 
functionals can be computed. Also the two-point function for pure quantum 
gravity can be obtained. This object is of central interest when questions 
about the scaling of the theory are to be addressed. 

On the technical side we have reviewed matrix models, which have de- 
veloped a status of interest on their own beyond dynamical triangulation in 
string theory and in condensed matter theory. 

Since there is at present no successful way to find a continuum theory of 
quantum gravity in higher dimensions than two, it is natural to ask, whether 
such a theory could be defined as the scaling limit of the corresponding higher 
dimensional discretized theory. It turns out that dynamical triangulation can 
be defined in three and four dimensions and that these theories have a phase 
transition similar to the two-dimensional case. While there has been some 
analytical progress in the understanding of entropy bounds on the number of 
triangulations in higher dimensions [86] most work in this field is numerical. 
Questions about the nature of the phase transition in four dimensions are 
not settled yet and exciting research remains to be done. 



Chapter 3 

The failure of quantum Regge 
calculus 



In the continuum path integral formulation of quantum gravity we are in- 
structed to compute the integral 



vol(Diff) 



(3.1) 



over diffeomorphism classes of metrics weighted with the exponential of the 
Einstein-Hilbert action, as has been discussed in chapter 1. Two different 
discretization schemes have been suggested. One, dynamical triangulation, 
provides a regularization of the functional integral (3.1). An explicit cutoff is 
introduced and one sums over all (abstract) triangulations with equilateral 
d-simplices of a d-dimensional manifold M, compare chapter 2. 

A different attempt to discretize the functional integral (3.1), called quan- 
tum Regge calculus (QRC), has been suggested, see [87, 88, 89] for reviews 
and references. This formalism is closely related to the classical coordinate 
independent Regge discretization of general relativity: One fixes a (suitably 
chosen) triangulation, while the dynamical degrees of freedom are the L link 
lengths. Thus the functional integration in (3.1) is replaced by 



f|dU J(l b ... ,10 6(A), 



(3.2) 



i=1 



where J(li,... ,l L ) is the Jacobian of the transformation dg^(£,) — > dl-t. 
The integral in (3.2) is over all link length assignments consistent with the 
triangle inequalities, as denoted by the delta function 6(A). While it provides 
a discretization of the integration in (3.1), this replacement does not provide 
a regularization of (3.1). Contrary to dynamical triangulation, no cutoff has 
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to be introduced [90]. One can of course choose to work with a cutoff which 
has to be taken to zero at the end of the computations. This will not make 
any difference in the argumentation below. 

The Jacobian in (3.2) is very complicated and its form is presently un- 
known. For analytical and numerical investigations it is usually replaced by 
local measures of the form Oi=i f(U)- 

If a sensible continuum limit of this theory existed, it would be ap- 
proached by taking the number L of links in the triangulation to infinity 
at a critical point of the statistical ensemble denned by (3.2) and by the dis- 
cretized Einstein Hilbert action. For this limit to make sense, the resulting 
continuum theory should not depend on details of the discretized theory like 
the chosen fixed triangulation or the local measure f(U). 

It is one central objective of this chapter to show analytically that such a 
continuum limit of quantum Regge calculus in its present formulation cannot 
be defined in any dimension d > 1 . 

Such a result has been indicated by simulations in two-dimensional quan- 
tum Regge calculus [91, 92, 93]. These simulations revealed that the KPZ- 
exponents (1.14) for the susceptibility of conformal matter coupled to con- 
tinuum quantum gravity could not be obtained. This has often been called 
the failure of quantum Regge calculus. Most discussion about this failure 
to reproduce continuum results in two dimensions has been centered around 
the choice of the measure. 

In the first part of this chapter some notation is introduced. In the 
second part we discuss the measures which have been suggested and used in 
the context of quantum Regge calculus. In the third part we show that in 
two dimensions all discussed measures for Regge calculus fail to reproduce 
continuum results. Not even the concept of length can be defined [2]. In the 
conclusion we discuss the situation in higher dimensions. We give a short 
review of an alternative approach, in which the continuum functional integral 
over metrics is restricted to piecewise linear metrics [94]. 

This chapter is partly based on work presented in [2]. 

3.1 Formulation of quantum Regge calculus 

Let us fix the connectivity of a triangulation with V vertices 1 , T triangles 
and L links of lengths li , . . . , l L of a two-dimensional closed manifold M 
with Euler-characteristic x- The number of links can also be expressed as 

1 We denote the volume of manifolds with the same symbol V. No confusion should arise. 
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L = 3V — 3x- We view the interior of the simplices as flat. As explained 
in section 2.1.1 the curvature of this piecewise linear space is concentrated 
at the vertices and each assignment of link lengths uniquely determines a 
metric on M. For a triangle this metric can be written as 

I xi +X2-X3) \ . , 

9j*v= ir , \ > ( 3 - 3 ) 

y {{^ +X2-X3) x 2 j 

where the squared edge lengths are defined as xt = I?, i = 1 , . . . , L. The 
area A of a single triangle can be expressed as 



d^v/gW = \ y/g = \ \j^x 2 -^[x A +x 2 -x 3 ) 2 . (3.4) 

Since in quantum Regge calculus the link lengths are the dynamical vari- 
ables one attempts to replace the functional integration over equivalence 
classes of metrics by an integration over all link length assignments which 
correspond to different geometries. Clearly, this replacement 



\X~V1 



Vol(Diff) 



J]dU J(l 1( ... , l L )6(A) = 



dM-(li,... ,k) (3.5) 



involves a highly non-trivial Jacobian J(li, ... ,l L ). The integration is over 
all link lengths compatible with the triangle inequalities. But not all assign- 
ments of link lengths define independent geometries as can be seen in the 
flat case. All vertices can be moved around in the plane, changing the link 
lengths without changing the flat geometry. The Jacobian has to be such 
that this additional invariance is divided out of the integral. 

Thus, we obtain an L = 3V — 3x-dimensional subspace of the infinite 
dimensional space of equivalence classes of metrics. Quantum Regge calculus 
replaces the functional integration in (1.1) by an integral over this finite 
dimensional subspace. It is hoped that in the limit L — > 00 expectation values 
of observables converge in some suitable way to their continuum values. 



3.2 Regge integration measures 

Presently, the form of the Jacobian in (3.5) is not known. However, by ap- 
pealing to universality one might assume that the choice of this measure is 
not very important. A large number of measures have been suggested and 
tried out in numerical experiments to study whether quantum Regge cal- 
culus in its present form agrees with the results of continuum calculations 



3.2. REGGE INTEGRATION MEASURES 



61 



and dynamical triangulations in two dimensions. The result of these inves- 
tigations is negative. In this section we will define the most important and 
common of these measures and study some of their properties. 



3.2.1 DeWitt-like measures 

One possible way to construct a measure for quantum Regge calculus in 
terms of the link lengths is to repeat the construction in section 1.2.1 with 
the constraint that the deformations of the metric are allowed by quantum 
Regge calculus. For an alternative discussion see [89] where the correspond- 
ing metric is called Lund-Regge metric, refering to an unpublished work by 
Lund and Regge. However, since the starting point of our discussion is the 
DeWitt metric we rather call the measure a DeWitt-like measure. Note how- 
ever that it is not the DeWitt measure restricted to piecewise linear metrics. 
This will be given in section 3.4.2. 

For a single triangle the variation of the metric (3.3) in terms of the 6x t 
is given by 



5g 



5xi ^(5xi + 6x 2 

^(Sxi + 5x2 — 6x 3 ) 6x2 



6x 3 ) 



(3.6) 



Using this and the DeWitt metric, the scalar product (5g,6g) T , which has 
been denned in section 1.2.1 can be computed for a single triangle with area 
A. It turns out that it simplifies considerably when we use the canonical 
value —2 for the parameter C in the DeWitt metric: 



(6g,6g) T = 



= - (2Sg^V p §g aP + C6g^g^g^6g aP ) 
= -2Adet(6g^y) det(gH, for C = -2. 

Using (3.3) and (3.6) leads to 



(3.7) 



(6g,6g>- 



[6x1,6x2,6x3] 



1 

16A 



1 


-1 


-1 " 




6x1 


1 


1 


-1 




6x2 


1 


-1 


1 




6x3 



(3.8) 



For a general two-dimensional triangulation this line element is given by the 
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sum of (3.8) over all triangles: 



= [6x 1t ... ,6x L ] M 



Sx L 



(3.9) 



Here M is an L x L matrix. My is the element corresponding to the links 
U and lj. If both U and lj belong to the same triangle t with area A t the 
corresponding non-diagonal entry is 



Mi 



(3.10) 



All other off-diagonal entries are zero. It follows that M has four nonvanish- 
ing off-diagonal entries in each row. Diagonal entries are given by 



1 1 

Mix = IT- + IT 



(3.11) 



where ti and t2 are the two triangles which contain U. Prom the i'th row of 
M the product (A tl A t2 )~ 1 can be factorized. Since each triangle has three 
sides, this amounts to factorizing rik=i fr° m ^ ne determinant of M: 



detM 



M< 3 



= n* 

k=1 

=: P(Ai,... ,A T )n A 



At + A 2 -A 2 -A 2 -A] -Ai 



(3.12) 



k=l 



P(Ai, . . . , A T ) is a polynomial in the areas of the triangles which vanishes, 
whenever two adjacent triangle areas vanish. It follows directly that P is a 
highly nonlocal function of the areas, since each monomial of P must contain 
at least half of the triangles. 

Therefore the DeWitt-like integration measure for quantum Regge calcu- 
lus is given by 

(MIt,... ,l L ) = const x V /P( ^ 1 ' ••• '^ llj^ij dlj 5(A). (3.13) 



n k =i A k 



On a first glance there are similarities with the continuum measure (1.55). 
Using g(£,) ~ A 2 we see that the power of the "local areas" in both measures 
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equals —3/2. However, reparametrization invariance fixes the continuum 
measure completely, while in the discretized version (3.8) might be multiplied 
by the diffeomorphism invariant area factor A 1_2|3/3 , leading to the measure 

/ H W T 

d|x(l 1; ... ,l L ) = const x V'"p T IT l j dl j 5(A) - ( 3 - 14 ) 

rik=i A k j=i 

This generalization stresses the fact that the power of the areas is not fixed 
in quantum Regge calculus by any obvious principle since the areas are dif- 
feomorphism invariant quantities with no direct local interpretation in the 
continuum. 

The measure (3.14) is highly nonlocal and thus not suited for numeri- 
cal simulations. Below we will show that it does not reproduce continuum 
physics. 

3.2.2 The DeWitt-like measure in other dimensions than two 

We can perform the derivation of the DeWitt-like measure in one and in 
higher dimensions. A one-dimensional piecewise linear manifold consists of 
L links of lengths U (xi = l 2 ),x = 1,... ,L, which are glued together at 
their ends. The canonical metric is thus given by = (xt). Therefore the 
DeWitt metric gives: 

C ^ 

(6g ) 6g) T =(1+-)^x i ? (5x i ) 2 . (3.15) 
i=i 

Thus the DeWitt-like measure in one dimension is 

l dl 

du(l] , . . . , l L ) = const x 1 T — p. (3.16) 

i =1 I? 

In this case the measure is local. The factor (1 + ^r) factorizes as in the 
continuum. Note that in d = 1 we cannot use the canonical value C = — 2 
for C since the DeWitt metric would be singular. 

In higher dimensions the situation is more complicated. For a d-simplex 
the measure is derived from the scalar product 

(6g,6g) T = SxM d 6x, (3.17) 

where Sx is the d(d 2 +1 1 -dimensional vector associated with link length defor- 
mations of the d-simplex. As in two dimensions, the complete norm for the 
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triangulation can be written as a straightforward superposition of the matri- 
ces M d . However, in dimensions larger than two, the entries of this matrix 
M.to tal depend on the xi's. In general, the measure is given by 

dn(l 1( . . . ,l L ) = const x y/det M*5 tal Y[ UdU 6(A), (3.18) 

i=i 

where 6(A) now stands for the generalized triangle inequalities. As an exam- 
ple we have computed the measure for the simplest closed three-dimensional 
manifold which consists of two tetrahedra glued together along the six links. 
The result is 

1 6 

d(i(l 1( ... ,l 6 ) = const x - J^fUdlt 6(A), (3.19) 

where V is the 3-volume of the manifold. 

Note that also in four dimensions, where the DeWitt metric is simply 
FUeMn^Kv d-9n-v, the determinant of M^ 1 does not evaluate to a constant, 
not even for the simplest 4-geometries. 

3.2.3 The DeWitt-like measures for special geometries 

We have not been able to obtain a general closed expression for the polyno- 
mial P. However, in a number of special cases it is possible to find explicit 
expressions for the measure (3.13) which can immediately be generalized for 
the measure (3.14) [2]. 

The easiest case of a closed piecewise linear manifold are two triangles of 
area A glued together along the three links. The measure is 

1 3 

du(l!,l 2 ,l3) = const x — TTljdlj 6(A). (3.20) 

A? U 

A tetrahedron has six links and four faces. The determinant of the 6x6 
matrix M can be computed. The measure is given by 

V 4 A- 6 

dn(l 1( ... ,l 6 ) = const x r 1 ' TTljdlj 6(A). (3.21) 

rii=iAiti 

Several other more complicated examples are possible. As a last example 
we display the measure for a hedgehog geometry which consists of K hat-like 
building blocks such as the one shown in figure 3.1. This block consists of 
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Figure 3.1: Building blocks for the hedgehog geometry which consists 
of K of these blocks glued together along the links l Qk and lb k . 



two triangles with areas A 2 k-i and A 2 k respectively, glued together along the 
links l2k-i and 1.2k- The hedgehog geometry is built by gluing these elements 
together along the links lb k and l ak+1 successively. After this gluing the 
vertices x and y have the order 2K, while the other vertices are of order two. 
Although this geometry might seem a bit artificial it plays an important role 
in the analysis of the phase transitions of two dimensional quantum gravity 
coupled to matter fields [69] . The matrix M for this geometry can be written 
in a block diagonal form. Its determinant can be computed, resulting in the 
measure 



d[i{h,... ,1 



3Kj 



const x 0£ ljdlj , Yl (Azic-i + A 2k ) 1/2 



nzN .3/2ii 
k=1 ^k k=1 

K K 



X 



k=1 



k=1 



(3.22) 



For an even number of elements this measure vanishes whenever the product 
over the even areas equals the product over the odd areas. For odd K, the 
measure is always positive. Clearly, this points to a severe sickness of the 
measure (3.13). 



3.2.4 Commonly used measures 

As has been mentioned above, the DeWitt-like measures (3.14) are not suited 
for numerical simulations. Also, we have demonstrated, that the measure in 
quantum Regge calculus is not fixed as in the continuum case. Therefore it 
is natural to search for other, local measures for which a sensible continuum 
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limit can be denned. The principle of universality should ensure that all 
reasonable choices of the measure result in the same continuum theory. 

For numerical simulations in two-dimensional quantum Regge calculus, a 
number of local measures have been suggested [87, 95, 96]. In all cases these 
are of the form 

i L 

d|i(ll , • • • , 10 = y-rj ^ J! l ^ dU 5(A) ' t 3 - 23 ) 

for some parameter a or (3, or derived from this 2 . These measures are arrived 
at by intuitively translating the continuum measure 0^11^ as 
Y\i = i UdU and by replacing the product EUgmS^) 11 by the Product of the 
"local" areas Ylk=i ^\f- The most popular choices are the scale invariant 
measure Yli=i t 1 an d the uniform measure rii=i d*d ~ Yl\=: UdU. 

Motivation for the first of these measures is often suggested by an analogy 
to the scale invariant continuum-measure 

^9n-v = n SW"^ 1 EI d 9^M> (3-24) 

xG M (i<-v 

which has been advocated for four-dimensional quantum gravity in [97, 98]. 
However, this measure is not diffeomorphism invariant, as has been stressed 
in section 1.2.1, although the opposite has been claimed and "proved" by 
wrong formal arguments [98, 99]. The DeWitt measure is the only diffeo- 
morphism invariant measure for continuum quantum gravity. Therefore the 
measure (3.24) cannot replace the DeWitt measure in continuum quantum 
gravity and any translations to discretized versions should be avoided. 

Also the DeWitt-like measure (3.13) cannot be used to motivate (3.23) 
in dimensions higher than one, since (3.13) is highly nonlocal while the mea- 
sures (3.23) are local. Note that also in four dimensions, where the DeWitt 
measure is simply 

nn^w. (3.25) 

xG M n<-v 

the discrete DeWitt-like measure is very complicated and highly nonlocal 
and does not equal the uniform measure (3.23) with <x = — 1 and (3 = 0, not 
even for the simplest 4-geometries. This clearly shows that translations like 

T L 

n gw a n d 9^w -> n A f ff n i ^ s(a) 

XGM h<-v j=1 i=1 

2 Sometimes the product over all triangle areas is replaced by the product over all areas 
assigned to the links. In our arguments this is not important and corresponds to choosing 
another value for the parameters a and |3. 
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or similar are far too simple minded, since the Jacobian of this variable 
transformation is not taken into account. The correct treatment of the dif- 
feomorphisms and the Jacobian leads to the results in section 3.4.2. 

However, it has been argued that, appealing to universality, the precise 
form of the measure need not to be known. In the continuum limit, av- 
erages of observables should result in their continuum values independent 
of details of the discretization like the exponents a and (3 in (3.23) or the 
underlying connectivity of the piecewise linear space. But numerical simula- 
tions revealed that averages of observables in the context of quantum Regge 
calculus 

• do not reproduce results of continuum quantum gravity. In two dimen- 
sions not even the KPZ exponents (1.14) could be obtained. 

• show a significant dependence on the measure. 

• show a significant dependence on the underlying triangulation. 
These are serious flaws in the Regge approach to quantum gravity. 

3.2.5 Numerical simulations in two-dimensional quantum Regge 



Quantum Regge calculus has repeatedly been used for numerical simulations 
in two, three and four dimensions. However, its status as a theory of quantum 
gravity remained unclear even in two dimensions. To clarify this situation, 
two groups have recently performed large scale numerical investigations in 
two dimensions where it is possible to compare with continuum results. 

Bock and Vink [91] have performed a Monte Carlo simulation of two- 
dimensional Regge calculus applied to pure gravity. They added the term 



to the action. Such curvature square terms have been introduced in the 
context of four-dimensional quantum gravity to overcome the unboundedness 
of the Einstein-Hilbert action from below, which is induced by conformal 
fluctuations [100]. In two dimensions this term is the only nontrivial part of 
the action. Its discretized version can be written as 



calculus 



(3 d 2 £,v/g TZ 2 



(3.27) 



Jm 




(3.28) 
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h 


(sphere) 


1 (torus) 


2 (bitorus) 


y in Liouville theory 


-0.5 


2 


4.5 


y in QRC from [91] 


r~s^ 5 .5 


^2.0 


r~s^ 5 .5 



Table 3.1: Values for the susceptibility exponent y for pure two- 
dimensional gravity at different topologies as computed in Liouville 
theory compared to the values from the Monte Carlo simulation in [91]. 
Bock and Vink do not give error bars for their values. However, the 
deviation from the Liouville values is significant. 



where 5 t is the deficit angle at the vertex i and A(i) is the area assigned to 
the vertex i by a barycentric division of the triangle areas around i. The 
addition of this term to the action allows a scaling analysis which gives the 
susceptibility exponent y. This exponent can for example be defined as an 
entropy exponent for the partition function 

Z(V) oc V 3 - y e AcV , (3.29) 

where the dominant exponential is nonuniversal and the universal part is 
given by the subleading power correction. Continuum calculations show, 
that for pure gravity y is given by 

7 = 2-5(1 -h), (3.30) 

compare (1.14). It has been shown in [71, 101], that the addition of a small 
curvature square term does not change the value of y in dynamical triangu- 
lation. 

To be conservative, Bock and Vink used the scale invariant measure (3.23) 
with <x = 1 and (3 = for their simulations. This enabled them to com- 
pare their results with previous calculations [96], where agreement between 
quantum Regge calculus and the formula (3.30) has been reported. They 
computed the susceptibility exponent for the topology of a sphere (h. = 0), a 
torus (h = 1) and a bitorus, that means a sphere with two handles (h. = 2). 
The results are compared with the continuum values in table 3.1. For spher- 
ical and bitoroidal topology they found that y in quantum Regge calculus 
differs significantly from Liouville theory. That means that quantum Regge 
calculus does not reproduce the fundamental KPZ-result (3.30). Note that 
the torus topology is not well suited to test the KPZ formula (1.14). For 
h. = 1 y is independent of the conformal charge of the matter coupled to 
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gravity. Furthermore y takes the same value 2 in the classical case without 
metrical fluctuations (c = — oo) and in the quantum case. 

The results of Bock and Vink have first been confirmed [102] and later 
been questioned [103, 104] by Holm and Janke. Using a very careful finite 
size scaling ansatz for pure quantum gravity on a fixed random triangulation 
of the sphere, they arrive at the value y = —10(2) for the susceptibility 
exponent [93], in the same setup. This is in clear disagreement with the 
prediction from Liouville theory and with the computations in [91, 96]. 

Using an alternative measure similar to (3.23) with <x = and (3 = \ it was 
found in [102] that a scaling can no longer be observed and one encounters 
numerical difficulties. The observables did not reach equilibrium values. A 
similar change in the scaling behaviour was observed in [91]. The value of 
y is changed when the scale invariant measure is replaced by nl=i i~^i> f° r 
C 7^ 0. Furthermore it has been demonstrated in [93], that the results of 
quantum Regge calculus depend on the chosen fixed triangulation. 



3.3 The appearance of spikes 

In this section we will show, that two-dimensional quantum Regge calculus 
does not reproduce results of continuum quantum gravity with any of the 
proposed measures mentioned above [2]. 

Let us recall that the partition function G(V, R) for closed two-dimensional 
universes with fixed volume V and with two marked points separated a 
geodesic distance R has the asymptotic behaviour (1.33) 

G(V,R) ~ V" 1/ V(vw) 3 ) foi-JL^oo. (3.31) 

Thus the expectation value of any power of the radius R can be calculated 
as 

dr r n G(V,r) ~ V n/4 . (3.32) 



We can show that (3.32) is not fulfilled in quantum Regge calculus. 

The physical reason is, that contrary to continuum quantum gravity in 
two dimensions, the volume, or equivalently the cosmological constant, does 
not set an intrinsic scale for the theory. It follows that the concept of length 
has no natural definition in this formalism and that every generic manifold 
degenerates into spikes. 
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3.3.1 General proof of appearance 

Theorem 1 For any value of a or (3 in (3.14) an d (3.23) there exists an 
n such that for any link I in a given arbitrary triangulation 



for any fixed value V of the spacetime volume. Thus the average radius 
(R), or some suitable power (R n ), does not exist [2j. 

Proof: In the situation of figure 3.2, which displays a small part of the piece- 



Figure 3.2: Parametrization of the link lengths and the triangle areas 
in a part of the piecewise linear manifold M. The vertex v is of order N. 
A spike can be formed by making the links li , . . . , In arbitrary large. 
When the area of the surface is fixed, the link lengths In+i , . . . , 
have to be small of order 1 /I] . 

wise linear manifold M, the vertex v has the arbitrary order N. The links 
connected to v are labelled as li, . . . ,In, while the links opposite to v are 
labelled as l N+ i , . . . , 1-2N- We want to analyze the Regge integral in that part 
of the phase space where the vertex v is pulled to infinity and forms a spike 
while the volume of M. is held bounded. In this situation the link lengths 
li , . . . , In are large, while the link lengths l N+ i , . . . , I2N have to be of the 
order . 

First we analyze the measure (3.23) with (3=0. Let A be a large num- 
ber. Then li can be integrated freely between A and infinity, while the 




(3.33) 
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integration over I2, . . . , In is restricted by triangle inequalities. The inte- 
gral over I2, . . . , In gives a factor l 1 , 1 N)k (In+i • • • Wi) 1 ^ -1 after symmetrizing 
over In+i, • • • , Un- An additional factor In+i . . . Un comes from the triangle 
inequalities for the adjacent triangles which contain the links In+i, • • • , 1-2N- 
Taking these factors together results in 



dli I, 



-Net 



dlN+1 . • • dl2N (In+i • • • 1-2n)" 



(3.34) 



A 



which exists if 



cc < 3 - (3.35) 



If that is fulfilled, the integral over l N+ i . . . I2N can be performed, which gives 



A 



-Net 



a: 



N 



dl, I] 



1-3N 



(3.36) 



A 



This means that (ly) = 00 for all n > 3N — 2, which proves the theorem 1 
for the measure (3.23) with (3=0.. 

To discuss the general form of (3.23) we have to parametrize the areas 
Ai , . . . , A 2 n- In our situation these are given as the products of a small and 
a large link length if v is pulled to infinity. Thus the product rL=i A, rP gives 
an extra factor (li . . . In) _(3 (In+i • • • 12n)~ 213 - After integrating over l 2 , . . . , In 
and after taking the triangle inequalities into account we get: 



A 



dl, Ij 



-N(«+P) 



2N — 1 -jr. 

dl N+1 . . . dl 2N (In+i • • • 12n)~ 1 , 



(3.37) 



which exists ifa + 2(3<3 — ^7. Then the integrals can be computed to give 



dli h 



1+N(|3-3) 



(3.38) 



Thus for all n > N(3 — (3) — 2 the expectation value of \\ is infinite. 

For the measure (3.14) the same analysis can be repeated if one notes, 

that from the nonlocal factor 



p(a; 



2(3/3 



. ,Af /3 ) the product flLi l f ^ 



can be split of. It turns out that this doesn't change the result and we 
get again (3.38). This completes the proof. We denote (l n ) = 00 as the 
appearance of spikes. 
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3.3.2 Spikes in special geometries 

Clearly, if a discretized theory shall make any sense, the results should not 
depend on details like the chosen triangulation. In fact however, the proof 
of theorem 1 shows, that (l n ) depends crucially on the order of the vertices, 
to which the link I is attached. To illustrate this dependence of quantum 
Regge calculus on the underlying triangulation further, we analyze some 
special geometries in the same way as in the proof to theorem 1. For some 
geometries even the partition function is ill defined for some measures, while 
for other measures or geometries (l n ) might be finite for some values of n. 

Spikes in a hexagonal geometry 



/ \/ \/ k \/ \/ \ 


\ / \ A 












/ \6k+6 . 


lc m 




/ \'a m 





Figure 3.3: Parametrization of a hexagonal geometry; k / m. All 
centers of the cells, which are marked by thicker lines form spikes. 

We analyze a regular triangulation with torus topology in which all ver- 
tices have the order 6, consisting of K hexagonal cells and parametrized as 
in figure 3.3. We want to study the measure (3.23) in that part of the con- 
figuration space, where all centers of cells form spikes while the total area of 
the surface is held bounded to prevent exponential damping from the action. 
That means that the 6K link lengths U, i = 1 , . . . 6K are very large while the 
3K link lengths l a]c , l bk , l Ck , k = 0, . . . , K — 1 are very small, of order l^ 1 . In 
cell number k one link, say l 61c+ i can be integrated freely from a large number 
A to infinity, while the integration over lek+2, • • • , Uk+6 is then constrained 
by the triangle inequalities. Integrating these out and symmetrizing over the 



3.3. THE APPEARANCE OF SPIKES 



73 



small links l Qk , l b]c , l C]c , k = 0, . . . , K - 1 gives 

poo k f A k 

ncr p) l6k+l n^wi b jcj^- a - 2p di ak di bk di Ck . ( 3 . 39 ) 
ja k=i j ° k=i 

Thus the integral over l aic ,l-b k >^c k only exists if <x + 2(3 < |. Furthermore 

(I?) = oo if n > 7 + 3cc. 

For the measure (3.14) we extract a factor rii=i V from the nonlocal 
factor. Performing the analysis gives the bound (3 < ^ on (3. (l^) is infinite 
if n > 4 — (3. These results are summarized in table 3.2. 

Spikes in a 12-3 geometry 

Similarly we can analyze the integration measures for quantum Regge calcu- 
lus for a triangulation of the torus where two thirds of the vertices have order 
3 and one third of the vertices have order 12, see figure 3.4. We consider that 



Figure 3.4: Illustration of the 12-3 geometry. All vertices of order 3 
form spikes such that the total area of the surface is held bounded from 
above, to prevent exponential damping from the action. 



part of the configuration space where all vertices of order 3 form spikes. The 
bounds on the exponents <x and (3 in (3.23) and (3.14) are depicted in table 
3.2. The resulting bounds are sharper because the order of the vertices which 
form spikes is lower which leaves less geometrical restriction. For a < — 1 
the the average of a link length cannot be defined with the measure (3.23) 
(independent of (3). The same holds true for the measure (3.14) if (3 > 0. 
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Spikes in degenerate geometries 



One can get sharper bounds for geometries with vertices of order two, like 
the hedgehog geometry introduced in section 3.2.3. As explained in that 
section, these geometries do occur in phases of simplicial quantum gravity 
in the framework of dynamical triangulation and it would be unnatural to 
forbid their occurrence as basically local parts of very large triangulations. 

As a first example we analyze the situation for the geometry shown in 



Figure 3.5: A degenerate geometry which has two vertices of high 
order while all other vertices are of order 4. Spikes can be formed in 
several ways. In the text we analyze the behaviour in that part of the 
phase space where every second vertex of order 4 forms a spike. 



figure 3.5 which has two vertices of high order while all other vertices are of 
order 4. This is slightly more regular than the hedgehog geometry associated 
with figure 3.1. In table 3.2 we present the results from our analysis of the 
situation, where every second vertex of order 4 forms a spike. While with 
measure (3.23) (I 2 ) is undefined for oc < — 1, the DeWitt-like measure (3.13) 
does not even admit the computation of the average (I). 

For the hedgehog geometry these measures themselves are illdefined, see 
table 3.2. Measure (3.23) is illdefined for oc = — 1 (the uniform measure), 
while for instance for the scale invariant measure (3.23) with oc = 1 and 
|3 = expectation values of I 2 and higher are infinite. 
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Geometry Measure Bounds Values of n 

s.t. (l n ) = oo 



hexagonal, 
see figure 3.3 


(3.23) 
(3.14) 


a + 2|3 < f, 


-|<a 


u > 7 + 3a 
n > 4- (3 


12-3, see 
figure 3.4 


(3.23) 
(3.14) 


a + 2|3 < f, 
(3 < 2 


-|<a 


TL > | + jOC 

n>1-§ 


degenerate, 
see figure 3.5 


(3.23) 
(3.14) 


a + 2|3 < f, 

P ^ io 


-2 < a 


n > 4 + 2a 
n>2-f 


hedgehog, see 
section 3.2.3 


(3.23) 
(3.14) 


a + 2|3 < 2, 
|3 <0 


-1 < a 


n > 1 + a 

n>-§ 



Table 3.2: Bounds for the exponents a and (3 for the measures (3.23) 
and (3.14) and values of n such that (l n ) = oo for various geometries. 



3.4 Concluding remarks 

3.4.1 The appearance of spikes in higher dimensions 

Although the derivation in the last section has been given only for the two- 
dimensional case, we expect, that one has to deal with the same problems 
in higher-dimensional quantum Regge calculus. The measures (3.23) have 
obvious generalizations in higher dimensions. An extra problem in higher 
dimensions is the unboundedness of the action which can be overcome by 
adding a cutoff involving for instance a curvature squared term. However, 
if the cutoff is taken to zero, we expect the spikes to reemerge. Actually, 
this has been observed in numerical simulations [105, 106]. In [106] four- 
dimensional quantum Regge calculus is explored on a fixed but not regular 
triangulation of the 4-torus. The triangulation is constructed from a regular 
triangulation with 31 vertices following [107] by inserting three additional 
vertices with low order. While in the conventional regular triangulation of 
the 4-torus all vertices are of order 30, the additional vertices are of order 5. 
The Einstein-Hilbert action plus a cosmological constant term is discretized 
in the usual way. A cutoff is introduced by the requirement that the local 
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fatness 

4> s = 576 (3.40) 

for every 4-simplex s is held bounded away from zero, 4> s > f = const > 0. 
Here we have introduced the 4- volume V s of the simplex s. The fatness is 
maximal for equilateral simplices and vanishes for collapsing ones. 
For the integration over the link lengths U the measure 

L 

d l i(l 1 ,...,l L )=n^" 1 dU5(A) (3.41) 

has been chosen for a between and 1 . In figure 1 of [106] the expec- 
tation value of I? is depicted as a function of the cutoff f at the values 
f = 2 m 10~ 6 , m = 9, . . . ,0 for a = 1, a = \ and a = ±. For a = 1 and 
ct = \ this expectation value tends to infinity if the cutoff is removed while 
for a = the situation cannot be decided from the numerical data in [106]. 
From the figure we deduce that the behaviour of (I?) is a power of f : 

(If) ~ r°- 2G[1 \ for o-= 1, (3.42) 
(if) ~ f-0-i6(D ) for o-= 1 (3.43) 

This clearly demonstrates the appearance of spikes in four dimensions and 
shows that one has to expect the dependence of the physical results on the 
choice of the measure and on the choice of the fixed triangulation in four 
dimensions as well as in two dimensions. 



3.4.2 The role of diffeomorphisms in two dimensions 

While we have shown that quantum Regge calculus in its present form is 
not a suitable candidate for a theory of quantum gravity it is an interesting 
question, whether one can develop a theory of quantum gravity in the spirit 
of quantum Regge calculus. 

A promising alternative approach to the functional integration (1.1) was 
studied in [108, 94, 109]. The integration over all metrics modulo diffeomor- 
phisms is replaced by the continuum functional integral over piecewise linear 
metrics with a finite number of singularities modulo diffeomorphisms. The 
singularities are the V vertices of the piecewise linear space at which the 
curvature is located. 

This approach is remarkable since at each stage in the calculations diffeo- 
morphisms are treated exactly. Thus it includes all metrics which are related 
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to a piecewise linear metric by diffeomorphisms. It turns out, that this is 
sufficient to reproduce a discretized form of the Liouville action. Therefore 
it is believed, that for V — > oo this discretized version approaches the contin- 
uum expressions in some sense. No ad hoc assumption about the form of the 
measure has to be made. Note that the connectivity of the piecewise linear 
space is not fixed. In particular, all metrics of dynamical triangulation are 
contained in this approach. 

We will give a short review of these ideas for the two-dimensional case. As 
in the continuum calculations, the starting point is the DeWitt metric. Since 
the reduction of the degrees of freedom involves only geometries but not dif- 
feomorphisms, the gauge fixing which leads to (1.72) can be performed in the 
same way as in the continuum. Now the integration V<$ over the conformal 
factor has to be restricted to those conformal factors which describe Regge 
geometries. For spherical topology, to which we confine the discussion here, 
there are no Teichmiiller parameters and one can choose a unique background 
metric §. Menotti et al. adopt the usual way of a stereographic projection 
of the piecewise linear surface on the plane with = 6^ [110]. Then the 
conformal factor for piecewise linear geometries can be parametrized as 

v 

e * = e 2Ao -Q| z _ z .|2(a t -1)_ ( 3 _ 44 ) 

Here 2noti is the angular aperture at the vertex i while the zt are the coordi- 
nates of the singularities in the complex plane. The sum of the deficit angles 
1 — <Xi has to equal the Euler characteristic, 

v 

2j1-a0=2, (3.45) 

that means that <xy can be expressed by the other angles. Ao is an overall scale 
factor. Note that this parametrization simply stems from a generalization 
of the common Schwarz-Christoffel mappings which map the complex upper 
half plane to an arbitrary polygon [111]. One advantage of this parametriza- 
tion over a parametrization via link lengths is that there are no triangle 
inequalities. 

Due to the presence of six conformal Killing vectors, compare section 
1.2.2, the gauge fixing is not complete and the conformal factor cf> is invariant 
under an SL(2,C) transformation given by 

, _ az + b ad _ bc = -| a.b.c.deC. (3.46) 
cz+d v ' 
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A short calculation gives 

(|>'(z\Ao,Zi,ai) = <J)(z',A£,z[,ai), (3.47) 

with 

v 

Ao = A -^(1 -ai)log|czi + d|. (3.48) 

i=1 

Geometrical invariants as the angles cxi and the area 

v 



A = e 2A ° 



d 2 z [ Jlz-ztl 2 *^- 1 ' (3.49) 



are invariant under the transformation (3.46). A counting shows, that as 
in the usual parametrization in Regge calculus there are 3V — 6 dynamical 
degrees of freedom. 

To compute the Fadeev-Popov determinant for these conformal factors 
one employs the same techniques as in the continuum. Namely one computes 
the change of 

det'(P+P) 

log - , — (3.50) 

det(cu Q , cu b )v 

under a variation of the parameters A , z t and at and integrates the result 
back. One gets [94]: 

(3.51) 

— -loglzi-Zjl + Ao^ai ) -^F(ai) 

X X 

This is the discretized form of the Liouville action. In the limit V — > oo 
it approaches the continuum Liouville action. The function F depends only 
on the angles ol\. For V — > oo, YjL=-\ F(<*i) goes to a topological invariant. 
The discretized Liouville action is invariant under the SL(2, C) transforma- 
tion (3.46). Indeed, it has been shown in [112] that this invariance induced 
by the conformal Killing vectors is sufficient to derive the form (3.51) of 
the Liouville action up to the sum over the F(<x0. Note that the discretized 
Liouville action is strongly unlocal. Overcounting of equivalent locally flat 
geometries is avoided by expression (3.51). Indeed, an angle oti = 1 simply 
does not contribute to (3.51). 
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The integration measure for the conformal factor is defined by the scalar 
product ds 2 = Jd 2 z e* ScjjScf). The change to the variables (qi, . . . , q3v) = 
((zi) x , . . . , (z v ) y , Ao, (Xi, . . . , <Xv-i) involves a Jacobian: 

V V-1 

^e*s,> = n d2zi Yl d(X i dAo V 7 ^!- (3-52) 
i=i j=i 



which is given by the determinant of 

h = 



with the conformal factor <b defined by (3.44). The diagonal entries Ju for 
i = 1 , . . . , 2V of this matrix require some regularization. The explicit form of 
the measure is unfortunately unknown even for the simplest geometries which 
makes the approach presently inaccessible for numerical tests. While it has 
been possible to extract the dependence of the Jacobian on the scale factor 
Ao, almost nothing is known about the other variables. It is an interesting 
open problem to show that the conformal anomaly arises from this measure. 



3.4.3 Conclusion 

We have shown that in two-dimensional quantum Regge calculus the expec- 
tation value of a suitable power of a single link length diverges. No natural 
scale is set by the measures (3.23) or (3.14). Even if the total area is kept 
bounded from above, the probability of having arbitrary large link lengths is 
not exponentially suppressed. In two-dimensional quantum gravity however, 
the probability of having two points separated a distance R is exponentially 
suppressed as e~ Rv/ ^ for a given cosmological constant A and suppressed as 
e -R 4 3 /v ^ area y Q £ spacetime is fixed. That means, that quantum 
Regge calculus in two dimensions does not provide a formulation of quantum 
gravity. In addition, the mere existence of the partition function depends 
crucially on the chosen measure as well as on the chosen fixed triangulation. 
Clearly, this is a severe problem: if quantum Regge calculus was a realis- 
tic theory, these microscopic details should be unimportant and not affect 
physical expectation values. 

Furthermore, we have shown that in four dimensions spikes have been 
observed in numerical studies of quantum Regge calculus. Numerical inves- 
tigations demonstrate that also in four dimensions the results do depend 
crucially on the underlying triangulation and on the measure. Thus our 
conclusion in four dimensions is the same as in two dimensions: Quantum 
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Regge calculus does not provide a theory of quantum gravity. A sensible 
continuum limit cannot be denned. Not even the concept of length can be 
defined unambiguously. 



Conclusion 



In this work we have analyzed the fractal structure of two-dimensional quan- 
tum gravity. A central result is that the spectral dimension equals two for all 
types of conformal matter with central charge smaller than one coupled to 
gravity. Other aspects of the fractal nature of the spacetime have been intro- 
duced and discussed such as the extrinsic and intrinsic Hausdorff dimension, 
the branching ratio into minimal bottleneck baby universes and the scaling 
of the two-point function. We have compared our results with the latest 
numerical data. Numerical simulations are performed in the framework of 
dynamical triangulation. We have discussed this method and illustrated its 
power by a computation of an exact expression of the two-point function in 
two-dimensional quantum gravity. Furthermore we have discussed quantum 
Regge calculus which has been suggested as an alternative discretization of 
quantum gravity. We have shown that this theory does not provide a suit- 
able discretization of quantum gravity since it disagrees with the continuum 
results. This point constitutes another central result of this work. 

There are several open questions which have to be addressed in future 
work. It would be very interesting to develop a method which would allow 
the solution of Liouville theory. More modestly, the computation of the 
geodesic distance or similar geometrical objects would be interesting. The 
derivation of the intrinsic Hausdorff dimension is not very well understood. 
An alternative derivation of this result would provide further insight into the 
theory. Related to this is the question, why the fractal dimension seems to 
equal four for unitary matter coupled to quantum gravity in two dimensions. 
Of course the final triumph would be an analogous analysis of the fractal 
properties of the quantum spacetime of four dimensional quantum gravity in 
the continuum approach and in the discretized approach. 

While this work settles the question about the viability of quantum Regge 
calculus in its present formulation it would be interesting to develop different 
schemes to discretize quantum gravity. We have discussed one possibility 
in chapter 3. It would be interesting to study further properties of this 
approach. 
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